Partial Recovery of Quantum Entanglement 

Runyao Duan, Yuan Feng, and Mingsheng Ying 



Abstract — Suppose Alice and Bob try to transform an en- 
tangled state shared between them into another one by local 
operations and classical communications. Then in general a 
certain amount of entanglement contained in the initial state will 
decrease in the process of transformation. However, an interesting 
phenomenon called partial entanglement recovery shows that it 
is possible to recover some amount of entanglement by adding 
another entangled state and transforming the two entangled 
states collectively. 

In this paper we are mainly concerned with the feasibility 
of partial entanglement recovery. The basic problem we address 
is whether a given state is useful in recovering entanglement 
lost in a specified transformation. In the case where the source 
and target states of the original transformation satisfy the strict 
majorization relation, a necessary and sufficient condition for 
partial entanglement recovery is obtained. For the general case 
we give two sufficient conditions. We also give an efficient 
algorithm for the feasibility of partial entanglement recovery in 
polynomial time. 

As applications, we establish some interesting connections 
between partial entanglement recovery and the generation of 
maximally entangled states, quantum catalysis, mutual catalysis, 
and multiple-copy entanglement transformation. 

Index Terms — Quantum entanglement, Entanglement trans- 
formation, Partial entanglement recovery, Entanglement cataly- 
sis, Majorization. 



I. Introduction 

QUANTUM entanglement is a valuable resource in quan- 
tum information processing. It can implement some 
intormation processing tasks that cannot be accomplished 
classically. As a consequence, entanglement has been widely 
used in quantum cryptography [1], quantum superdense coding 
[2], and quantum teleportation [3]; see [4], Chapter 12 for an 
excellent exposition. Due to the great importance of quantum 
entanglement, a fruitful branch of quantum information the- 
ory named quantum entanglement theory is currently being 
developed. 

Since quantum entanglement exists between different sub- 
systems of a composite system shared by spatially separated 
parties, a natural constraint on the manipulation of entangle- 
ment is that the separated parties are only allowed to perform 
local quantum operations on their own subsystems and to 
communicate to each other classically (LOCC). Using this 
restricted set of transformations, the parties are often required 
to optimally manipulate the entangled state. One of the central 
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problems of quantum entanglement theory is thus to find the 
conditions for when an entangled state can be transformed into 
another one using LOCC. This problem can be approached in 
two different, but complementary, contexts: the finite regime 
and the asymptotic regime. In the asymptotic regime Bennett 
and his collaborators [5] proposed a reversible protocol which 
shows that any two bipartite entangled pure states with infinite 
copies can be converted into each other without any loss of 
entropy of entanglement. Since in practice one can only have 
finitely many copies of an entangled state, it is of great interest 
to consider the problem of entanglement transformation in a 
finite (non-asymptotic) setting. Arguably, the most important 
step in the finite regime was made by Nielsen in [6], where 
he reported a necessary and sufficient condition for a bipartite 
entangled pure state to be transformed into another pure one 
deterministically using LOCC. Suppose two distantly located 
parties, Alice and Bob, share an entangled state in C n ®C n 
with Schmidt decomposition \ip) = Yl7=i \f&iVA)\iB), where 
c*i > 0.2 > • • • > a n > are Schmidt coefficients and 
ai ~ \^ A ) anc ^ V B ) are orthonormal bases of Alice's 
and Bob's systems, respectively. Suppose the two parties 
want to transform \if>) into another state \cp) with Schmidt 
decomposition \ip) = £)"=i a/5Ka)I*b)' where fa > j3 2 > 
■ ■ ■ > Pn > and Y2?=i P' 1 = The orthonormal bases 
\%a) and \i' A ) (similarly, and \i' B )) are not necessarily the 
same. For the sake of convenience, we write if> = (ai ,-•■,«„) 
and ip — (f3i, ■ ■ ■ , (3 n ) for the respective ordered Schmidt 
coefficient vectors of \ijj) and \<p). Nielsen proved that Alice 
and Bob can achieve this transformation of \if>) to \ip) with 
certainty using LOCC, written \ip) — > \ip), if and only if 
ip -< f- Here the symbol '-<' stands for 'majorization relation', 
and ip -< ip holds if and only if 

i i 
^ &i < ^ A for all 1 < I < n 

i=l i=l 

and Y^i=i a i = SiLi &i< If all inequalities in the above 
equation hold strictly and Y^h=x a i = SILi A> l h en we say 
that if) is strictly majorized by <p. Majorization is an interesting 
and well-developed topic in linear algebra. For more details, 
we refer the reader to [7] and [8]. 

Nielsen's theorem establishes a connection between the 
theory of majorization and entanglement transformation. It is 
of fundamental importance in studying entanglement transfor- 
mation and has many interesting corollaries. For example, by 
taking limits the asymptotic result of Bennett et al can be 
recovered from Nielsen's theorem. Unlike the transformations 
in the asymptotic regime, a direct implication of Nielsen's the- 
orem is that the amount of entanglement decreases during the 
deterministic transformation. Let E(\if>)) = — Y17=i ai 1°§2 a i 
be the entropy of entanglement of \ifj) , Then by the properties 
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of majorization, it follows that \ip) — > |<^) implies _E(|</?)) < 
-E(IV')) [6]. Indeed, these properties of majorization imply 
that any well-behaved entanglement measures, such as Renyi's 
entropy, or any other suitable concave functions, also decrease 
under LOCC. Intuitively, this means that a certain amount 
of entanglement will be lost in a LOCC transformation. It 
would be desirable to save some entanglement lost and reduce 
the net loss of entanglement in the transformation, since the 
saved entanglement can be used, for example, to increase the 
classical capacity of a quantum channel [9]. 

The possibility of recovering lost entanglement was first ob- 
served by Morikoshi [10]. We outline Morikoshi's recovering 
scheme as follows. Suppose Alice and Bob share an entangled 
state and they can transform it into \<p) by LOCC. As we 
mentioned above, this process is generally a dissipative one in 
the sense that the quantity of entanglement in the target state 
is less than that in the source state. Suppose now an auxiliary 
state |x) is supplied to Alice and Bob. Instead of transforming 
\ip) into \ip) directly, they perform collective operations on 
the joint state \ip) (g> |x), and transform it into another joint 
state \<p) ® \ui). Of course, as required by Nielsen's theorem, 
entropy of entanglement of the whole system decreases too. 
But by choosing a suitable auxiliary state |x), sometimes a 
state \u) with more entropy of entanglement can be obtained. 
Intuitively, this process enables part of entanglement lost in the 
original transformation to be transferred to the auxiliary state, 
and it was termed partial entanglement recovery. Morikoshi 
demonstrated that partial entanglement recovery for a trans- 
formation between 2x2 states is always possible by using a 
2x2 auxiliary state. 

Partial entanglement recovery for transformations between 
higher dimensional states was considered by Bandyopadhyay 
et al. in [12]. To avoid trivial cases (a perfect recovery can 
always be achieved by letting |x) = \<p) and \u>) = \ip)), a 
notion of genuine partial recovery was introduced. A partial 
recovery scheme is genuine if the dimension of the auxiliary 
state is smaller than that of the original source state. Then it 
was proven that for any states \ip) and \ip) such that ip is strictly 
majorized by <p and n > 2, a genuine partial recovery is always 
possible by using only 2 x 2-dimensional auxiliary states. This 
extensively generalized the result in [10]. The possibility of 
genuine partial recovery for the transformation of \ip) to \<p) 
such that ip is not strictly majorized by ip was also examined 
carefully in [12]. However, several fundamental problems 
concerning partial entanglement recovery are still open. For 
example, the existence of genuine partial recovery in the case 
of a n — fi n is unknown. Furthermore, the proof confirming the 
existence of auxiliary states for partial entanglement recovery 
presented in [12] is nonconstructive. In general, this proof 
method does not provide a way to find these auxiliary states 
efficiently. 

In this paper, we study the feasibility of partial entanglement 
recovery. We consider the problem of whether a given entan- 
gled state can be used to recover some entanglement lost in 
a specified transformation. Our motivations are twofold. The 
first one is more theoretical. In some sense, the process of 
partial entanglement recovery reveals a new kind of application 
of quantum entanglement: it can be used to store some entan- 



glement lost in information processing tasks. So it is of great 
theoretical interest to characterize the entanglement recovering 
ability of a given entangled state, as it may lead to a better 
understanding of some fundamental properties of quantum 
entanglement. In addition, as will see later, the solution of the 
above problem leads us to a rich mathematical structure and 
provides new insight into the process of partial entanglement 
recovery. The second motivation is more practical. Suppose 
we are required to perform a couple of different entanglement 
transformations. In most applications the available entangled 
states shared between two parties are pre-specified and very 
limited. A solution to the above problem will help us to 
determine whether partial entanglement recovery for these 
transformations is possible with other pre-specified entangled 
states. It is also worth noting that this problem is more general 
than the ones discussed in [10] and [12], and its solution 
automatically resolves many trivial cases. 

To state the above problem more formally, let us assume 
that \ip) and \<p) are the source state and the target state of 
the specified transformation, respectively, and let |x) be the 
given auxiliary state. Furthermore, suppose that \ip) can be 
transformed into \ip) with certainty using LOCC. Our goal is to 
determine whether there exists another state \u>) satisfying (i) 
the transformation of |^)<S>|x) to l</?)® I w ) can be implemented 
with certainty using LOCC, and (ii) \u>) is more entangled 
than |x). Next, we clarify a subtle point: what is meant by 
the statement that a state is more entangled than another one? 
An exact mathematical definition is needed. One way to do 
this is to use some measures of entanglement such as entropy 
of entanglement mentioned above, as in [10] and [12]. Note 
that for deterministic transformations, a single measure of 
entanglement is usually not enough to quantify entanglement 
amount since there exist incomparable states \ip) and \ip), i.e., 
neither \ip) — > \ip) nor \ip) — > \ip) is possible [6]. So in the 
present paper, we adopt an alternative view-point: we say that 
\ui) is more entangled than |x) if \u) — > |x) and |x) \oj). By 
Nielsen's theorem, this is equivalent to u ~< \ an d X 00 (here 
both x and ui are in non-increasing order). We believe that this 
view-point is more reasonable than only considering a single 
measure. Now the mathematical problem of the feasibility of 
partial entanglement recovery can be clearly formulated as 
follows: 

Problem 1. Given a triple of states (\ip), \ip), |x)) such that 
ip -< tp, determine whether there exists a state \u) such that 

ip^X^f®*^, u)<x> an d X UJ - 

If such a state \u>) does exist, then we call it a solution of 
Problem 1. In the above formulation we made no additional 
assumptions on the dimension of \ except that it is finite. 
So even in the case that the dimension of x is larger than or 
equal to that of ip (and tp), the above problem still makes sense. 
This enables us to consider the process of partial entanglement 
recovery in a general mathematical framework. It is also worth 
pointing out that Problem 1 cannot be directly solved by linear 
programming methods because the majorization relation ip ® 
X -< <p ® oj cannot be expressed by linear constraints, unless 
we know how to order a tensor product of two probability 
vectors, tp®u. The main difficulty here is that the order of 
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cp ® u/ is not related in any simple way to the orders of cp and 

10. 

The principal aim of the present paper is to solve Problem 1 
stated above. We first introduce three indices of uniformity for 
bipartite entangled pure states. With the aid of these indices, 
we prove that whether \x) can save some entanglement lost 
for the transformation of l^) to \cp) only depends on the target 
state \cp) and the presence of the equalities in the majorization 
relation ip -< cp. To be concise, let be a state with m 
distinct Schmidt coefficients, say, > ■ ■ ■ > a' m > 0. If 
m > 1, then the maximal local uniformity of denoted by 
L u {\ip)), is given by the maximal ratio of a' i+1 and a[ for 
all 1 < i < m — 1. In contrast, the global uniformity of \ip), 
denoted by g u (\ ip)), is given by the ratio of a' m and a[. In 
the special case of m = 1, both indices are defined to be 1. 
These indices have many useful properties. Indeed, they are 
key tools in studying partial entanglement recovery. With these 
notions, Problem 1 is completely solved in the case where ip is 
strictly majorized by cp (Theorem |3.21 >. We achieve this goal by 
considering two cases. First, Problem 1 is examined carefully 
for a special case where all nonzero Schmidt coefficients of 
\x) are identical, i.e., L u (\x)) = or L u (\x)) = 1- Second 
we consider the general case that < L u (\x)) < 1 and prove: 

(1) if L u (\x)) > gu{\f)), then \x) can recover some 
entanglement lost for the transformation of \ip) to \cp); 

(2) if L u {\x)) = g u (\cp)), then there is only a special form 
of \cp) for which |x) can save some entanglement lost in the 
transformation of 1^) to \cp); and 

(3) if L u (\x)) < g u (\cp)), then |x) cannot recover entangle- 
ment lost in any transformation with the target \cp). 

It should be pointed out that the proof we present provides 
an explicit construction of the resulting state \uj). In view of 
this, the above results are very useful in pursuing practical 
applications of partial entanglement recovery. Some interesting 
special cases of these results are also discussed. 

For the case where ip is not strictly majorized by cp, a com- 
plete solution of Problem 1 appears to be very difficult. Nev- 
ertheless, two sufficient conditions for partial entanglement 
recovery are presented (Theorems 14.11 and 14.21 1. Employing 
these conditions as tools we show that the genuine partial 
recovery is not always possible when the dimension of the 
target state is larger than 2x2. For example, if n = 3 and 
ct\ = (3i, \x) should be at least a 3 x 3 entangled state, which 
means any recovery scheme cannot be genuine. (This result in 
fact has been obtained implicitly in [12]). When a\ = f3\ and 
an = fin, we show that 4 x 4-dimensional auxiliary states are 
necessary and sufficient. In the case where n = 4, a genuine 
partial recovery is not possible. On the other hand, even in 
these special cases, it still makes sense to consider whether |x) 
is useful in recovering entanglement lost in the transformation 
of |^) to \cp). 

Besides the mathematical characterization of partial entan- 
glement recovery outlined above, we also present an algorith- 
mic approach to Problem 1 . Let n and k be the dimensions of 
ip (as well as cp) and x ( as we U as w )> respectively. Our goal 
now is to design polynomial time algorithms in n or/and k 
to solve Problem 1. As mentioned above, the main difficulty 
in solving Problem 1 lies in the fact that the order of the 



tensor product cp ® to cannot be determined by a simple 
method even after we know the orders of cp and to. Thus 
one cannot apply standard linear programming techniques 
directly. A naive enumeration of the possible orders of <p ® to 
yields about (nk)\ results, which is intractable. A simple but 
powerful lemma is introduced to reduce the number of orders 
of the tensor product. The basic idea behind this lemma comes 
from the observation that for a fixed cp, ip ® u> has at most 
0((fcn) 2 ( fc-1 )) different orders when ui varies. This number of 
the possible orders is only a polynomial in n when k is treated 
as a constant. For each possible order, we can employ linear 
programming methods to solve the majorization inequality 
ip<8>X ~< cp&tu. Consequently, an algorithm of time complexity 
0(n 2fc_1 log 2 n) is obtained (Theorem l5.11 l. This algorithm is 
not efficient in the case where k can vary freely. Fortunately, 
by examining the mathematical structure of partial entangle- 
ment recovery carefully, we can further refine this algorithm 
into a new one with time complexity 0(n 2 k 4 ) (Theorem 15. 2\ . 
Therefore we can efficiently determine the feasibility of partial 
entanglement recovery by using algorithmic methods. 

To illustrate the utility of the above results, we show 
that partial entanglement recovery also happens in situations 
such as quantum catalysis, mutual catalysis, and multiple- 
copy transformation. As an interesting application, we con- 
sider the generation of maximally entangled states using the 
scheme of partial entanglement recovery. We prove that any 
transformation with the Schmidt coefficient vector of the 
source state being strictly majorized by that of the target state 
can always concentrate some partially entangled state into a 
maximally entangled one. We also find a close connection 
between partial entanglement recovery and quantum catalysis 
(see [18], [19], [21]). That is, if a transformation can be 
implemented with certainty by using some quantum catalyst, 
then entanglement lost in the transformation can be partially 
recovered by a suitable auxiliary state. Moreover, we show that 
partial entanglement recovery is directly connected to mutual 
catalysis [23]. As a consequence, a systematic construction of 
the instances with mutual catalysis effect is sketched. When 
we consider the possibility of partial entanglement recovery 
in multiple-copy transformations (see [15], [19], and [21]), we 
notice a very interesting phenomenon: although an auxiliary 
state cannot be used to do partial entanglement recovery for a 
single-copy transformation, it can recover some entanglement 
lost in certain multiple-copy transformations. 

The rest of the paper is organized as follows. Section II 
presents some notations and concepts, including the definitions 
of uniformity indices. In Section III, we present a complete 
solution to Problem 1 in the case that ip is strictly majorized by 
cp. We consider general transformations in Section IV and give 
two sufficient conditions for partial entanglement recovery. 
Some special but interesting cases of these conditions are 
investigated in detail. In Section V we discuss the feasibility of 
partial entanglement recovery from an algorithmic viewpoint 
and present two algorithms to solve Problem 1 . To understand 
whether partial entanglement recovery is possible in situations 
such as quantum catalysis, mutual catalysis, and multiple- 
copy transformation, we give more examples and discussions 
in Section VI. In Section VII, we draw a brief conclusion. 
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The proofs of some lemmas and theorems are completed in 
Appendices. 

II. Preliminaries 

First, it is helpful to introduce some notations associated 
with finite dimensional vectors. Let x be an n-dimensional 
vector. The dimension of x is denoted by dim(x), i.e., 
dim(x) = n. The notation x* is used to stand for the vector 
that is obtained by rearranging the components of x into 
non-increasing order. Similarly, x 1 ' denotes the vector that 
is obtained by rearranging the components of x into non- 
decreasing order. The notation x® k denotes the direct sum of 
x with itself k many times. In particular, for constant c, c® fc 
is the fc-dimensional vector (c, c, . . . , c). If every component 
of x is nonnegative, then we can write 

X (Xi , . . . , X m J, 

where x\ > • • • > x' m > 0, fc, > 1 for each i = 1, . . . , m, 
and X)2=i = n - The above form of x* is usually called the 
compact form of x. It is obvious that the compact form of a 
nonnegative vector is unique when the dimension of the vector 
space under consideration is fixed. 

The sum of the to largest components of the vector x is 
denoted by e m (x). That is, e m {x) — YHL\ X \- ^ i s easily 
to verify that e m (x) is a continuous function of x for each 
m = 1, . . . , n. 

We say that x is majorized by y, denoted by x -< y, if 

e m (i) < e m {y), for all 1 < m < n — 1 (1) 

and e n (x) = e n (y). If all inequalities in Eq. ([TJ are strict and 
£n{x) = e n (y), then we follow the terminology in [12] and 
say that x is strictly majorized by y, denoted by x < y. 

A vector x is a segment of a vector y if there exist i > 1 
and k > such that x = (y i: y i+1 , y i+ k)- 

Now we apply the above terminology to bipartite entangled 
pure states. Let \ip) be an n x n entangled pure state with 
ordered Schmidt coefficients otx > a% > • ■ ■ > a n > 0. 
As we have mentioned in the introduction, the symbol ip is 
used to denote the Schmidt coefficient vector of \ip), i.e., 
ip = (a%, . . . , a n ), which is just an n-dimensional probability 
vector. We often identify the compact form of ip with the 
compact form of state \ip). We call \ip) an n x n maximally 
entangled state if the compact form of \ip) reduces to (1)©"; 
otherwise we say that \ip) is a partially entangled state. If tp'^- 
is a segment of tp*, then we call \tp') an unnormalized state. 

To apply Nielsen's theorem to unnormalized states, it can 
be restated as: \tp) — > \<p) if and only if ip -< tp. 

We define S(\<p)) to be the set of all n x n entangled 
pure states \ijj) which can be directly transformed into \<p) 
by LOCC. By Nielsen's theorem, S(\tp}) = {\ip) : ip -< tp}. 
We also define S°(\ip}) to be the set of all n x n states 
\ip) such that ip is strictly majorized by tp, i.e., S°(\tp)) = 
{\ip) : ip < tp}. It should be noted that states \ip) in S(\tp)) 
are required to have the same dimension as \tp). Such a 
requirement forces us to distinguish S(\tp)) from S(\tp')) when 
\tp) and \tp') are essentially the same state but their dimensions 
are different. For example, let tp = (0.5, 0.25, 0.25) and 



tp' = (0.5,0.25,0.25,0). It is obvious that the states \tp) and 
\tp') are essentially the same. However, according to the above 
definitions, S(\tp)) is completely different from S(\tp')). This 
design decision in defining S(\tp)) enables us to considerably 
simplify the presentation of our main results. The same remark 
also applies to the definition of S°(\ tp)). 

In this paper, the phrase 'bipartite entangled pure state' is 
used frequently. So, for convenience, sometimes we abbreviate 
it to 'state' or 'quantum state'. This should not cause any 
confusion from the context. 

Now we introduce three notions which are key mathematical 
tools in describing partial entanglement recovery. 

Definition 2.1: Let \ip) be an n x n partially entangled state 
with compact form (a / 1 ® fcl , . . . ,aj® fcm ), where n = hi 
and m > 1. Then 

(i) the minimal local uniformity of \ip) is defined by 

iu(IV')) = min{^±i : 1 < i < to}; 

(ii) the maximal local uniformity of \ip) is defined by 

L u (\ip}) = max{^yi : 1 < i < m}; 

(iii) the global uniformity of \ip) is defined by 

It is easy to see that the minimal local uniformity, the 
maximal local uniformity, and the global uniformity of a 
quantum state \ip) with ip^ = (ai, . . . , a n ) may be rewritten 
in a slightly different way: 

l u (\iP)) = min{^±i : 1 < i < n}; 

L u {\ip)) — max{ 1+1 : 1 < i < n and ai > a^+i}; 
at 



The above rewriting will help us to simplify some proofs. 

From the above rewriting of Definition 12. U it is easy to see 
that both l u (\ip)) and g u (\ip)) are continuous with respect to 
\ip). Thus it is reasonable to define the minimal local unifor- 
mity and the global uniformity of a maximally entangled state 
to be 1. However, such a continuous property does not hold 
for the maximal local uniformity. To keep many properties 
of these indices valid even in the case that the quantum state 
under consideration is maximally entangled, it is convenient to 
define the maximal local uniformity of a maximally entangled 
state to be 1. Also, for the sake of convenience, when the 
dimension of the state under consideration is one-dimensional, 
we define the uniform indices as 1, 

In applying the above definitions of uniformity indices, 
it should be noted that the dimension of \ip) is somewhat 
arbitrary, as one can append zeroes to the vector ip and thereby 
increase its dimension without changing the underlying quan- 
tum state. Suppose that the number of nonzero components of 
ip is n. If \ip) is treated as an n x n state, all the above three 
uniformity indices are positive. However, if we append zeroes 
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to ip and yield a state \ip'}, then the uniformity indices of \tp') 
are changed rapidly. For example, let ip = (0.5, 0.25, 0.25) 
and ip' = (0.5,0.25,0.25,0). It is obvious that both the 
minimal local uniformity and the global uniformity of \ip) are 
0.5. However, the minimal local uniformity and the global 
uniformity of \ip') are changed into 0. To avoid any confusion 
that may be caused by the phenomenon that we just mentioned 
in the above definition, the dimension of the states are treated 
as fixed. In other words, if \ip') is obtained from \ip) by 
appending zeros in its Schmidt coefficient vector, they may 
be thought of being two different states. Therefore, it is 
reasonable to allow that sometimes l u (\ip)) ^ l u (\ip')) (as 
well as L U (|V» ? L u (\f)) and g u (\iP)) £ g u (\ip'))). 

Some simple but useful properties of the three indices 
defined above are presented in the following: 

Lemma 2.1: Let \ip) and \<p) be two quantum states 
with compact forms ip^ — (a'® kl , . . . , a'® kr ) and ip^ = 
(f3'® h , . . . , 0'® l ° ), respectively. Then 

(1) 0<W)),L u (\iP)),g u (\i>)) < 1. 

(2) v-Hm^guim^um- 

(3) gM) < L*- 1 ^)). 

(4) g u m)<W))<L a (\ip)). 

(5) if r = s and ol i = j3\ for i = l,...,r, then 

= k(\<p)). Similarly, L u (\ip)) = L u {\<p)) and 
9u(\ip)) = 9u{W))- 

(6) if \i>) -f \<p), then 9u (\tJ))) > g u (\<p)). 

Proof. (l)-(5) follow immediately from Definition I2.ll 
(6) follows directly from Definition I2.1I and the fact that if 
|V) -» \<p) then a[ < /3[ and a' r > f3' s . □ 

We give some remarks on the above properties. (1) shows 
that the three indices of minimal local uniformity, maximal 
local uniformity and global uniformity are all between and 1 . 
Moreover, they take value 1 if the state is maximally entangled. 
The minimal local uniformity and the global uniformity take 
the value if the state in question has zero as one Schmidt 
coefficient, while the maximal local uniformity takes value 
if it is a maximally entangled state in a state space with lower 
dimension, i.e., with a compact form ((^)® m , 0®"~ m ) for 
some m < n. If l u (]ip)) = or l u (\ip)) = 1, i.e., \ip) has 
zero as a Schmidt coefficient or it is maximally entangled, 
then both the inequalities in (2) hold with equalities. In the 
case that < ^ U (|V)) < 1, the first equality in (2) holds 
if the distinct Schmidt coefficients of \ip) form a geometric 
sequence; while the second equality holds if \ip) has at most 
two distinct Schmidt coefficients. The equality in (3) holds 
if and only if the distinct Schmidt coefficients of \ip) form a 
geometric sequence. (4) can be analyzed similarly. (5) means 
that these indices only depend on distinct Schmidt coefficients 
of the state. (6) indicates that global uniformity is decreasing 
under LOCC. 

In addition to these trivial properties displayed in Lemma 
12.11 the following lemma presents three more interesting 
properties of global uniformity and minimal local uniformity: 

Lemma 2.2: Let \ip) and \ip) be two quantum states. Then 

(1) 9u(\ip) <S> \<f)) = gu(\ip))9u{\f))- In particular, 
g u (\i>)® m ) = 9u(\i>)) my m> 1. 

(2) Z„(|V> ® M) > min{Z„(|V>)), l u (\<p))}. 



(3) l u (m < lu(\i>)® k ) < min{^,^}, where ^ = 

Proof. (1) follows immediately by Definition 12.11 (3) is a 
simple application of (2) and Definition 12.11 So it is enough 
to prove (2). 

Let V X = (ai, • • • , a m ) and ip 1 = {[3 X , ■■■ , (3 n ). Since 
the minimal uniformity /„ is a continuous functional, we can 
assume without loss of generality that all components of ip 
and ip are positive. Let 

a = a p j3 q and b = a r (3 s , a < fe, 

be any two successive elements of the ordered probability 
vector (ip <g) ip)K It is obvious that r < m or s < n. Suppose 
r < rn is the case, let us try to prove 

urn < - b - (2) 

Indeed, from the definition of l u we have 

= <*r+f* < 1. V> 
OtrPs — 

Thus a r+ iP s < b. But, since a and b are successive elements, 
a r +i/3 s cannot belong to the interval (a, b), that is, 

a r+1 f3 s < a. (4) 

From Eqs. (0) and we get immediately Eq. (0. 

If r = to then we can be sure s < n. Thus we can apply 
analogous arguments to prove that 



In both cases one has 



uim{l u (\iP))MW))} < 

Since this is true for any successive a < b, we have proved 
statement (2). □ 

The above lemma deserves some more remarks. Intuitively, 
(1) shows that the global uniformity is multiplicative under 
tensor product. (2) means that the tensor product of two states 
is at least as uniform as one of them. (3) provides an upper 
bound and a lower bound respectively for the minimal local 
uniformity of any state consisting of multiple copies of a 
given state. More interestingly, it shows that the minimal local 
uniformity ofa2x2or3x3 state remains invariant under 
tensor products involving multiple copies. 

One of the most interesting applications of the uniformity 
indices introduced above is that they provide a characterization 
of when a strict majorization relation holds. 

Lemma 2.3: Let \ip) and |x) be two quantum states, and let 
S°(\(p)) <S)\x) denote the set of all states of the form \ip) ® |x) 
with |V) in S°{\ V )), i.e., S°(\cp)) ® \ X ) = M ® |x> : W e 
Then 



S°(\ip)) ® |X) C S°(\<p) 
Proof. See Appendix A. 



\ X ))^lu(\x))>9u(\v)). 



□ 



Roughly speaking, the above lemma shows that if the aux- 
iliary state |x) is uniform enough, then the strict majorization 
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relation ip ® \ <1 V ® X can be kept providing ip<(p, and vice 
versa. What we would like to emphasize here is that the only 
constraint on the source state \tjj) is tp <] cp. 

In the introduction we have frequently used the notion of 
partial entanglement recovery. We present a rigorous definition 
as follows. 

Definition 2.2: Let \ip) and \ip) be two n x n states, and let 
|x) be a kx k state. We say that |x) can do partial entanglement 
recovery for the transformation of \t/j) to \tp) if there exists a 
kxk state |w) such that 

(i) both the transformations of \ip) ® |x) to |(/?) (8) |cj) and 
|w) to |x) can be realized with certainty under LOCC. That 
is, \ip) <g> |x) -> \<p) <8 |w) and |w) -> |x); 

(ii) the transformation of |x) to |w) cannot be achieved 
with certainty under LOCC. That is, |x) 

Some remarks follow: 

1) In the above definition, both the dimensions of the 
source state \ip) and the target state \<p) are n x n. 
Similarly, the dimensions of the auxiliary state |x) and 
the resulting state \u>) are both kxk. These constraints 
are reasonable since during the transformation process 
the state space under consideration isn't modified. In- 
tuitively, and \<p) are two different states of the 
same two particles entangled between Alice and Bob. 
The dimensions of these particles are assumed to be 
finite and fixed. Similar arguments apply to the states 
|x) and \oj). This is in fact the reason that we have 
to require that the dimension \ip) is fixed in defining 
S(\<p)), S°(\cp}), l u (\<p)), L u (\ip)), and since all 
of them are introduced in this paper to describe partial 
entanglement recovery. 

2) According to Nielsen's theorem, the above definition can 
be rewritten as: a k x k auxiliary state |x) can do partial 
entanglement recovery for a transformation of \ip) to \<p) 
if there exists another kxk state \ui) such that all three 
relations ip (g> x ^ <P ® w -< x> an d X 1 ^ ^ h°kl 
simultaneously. 

3) It is obvious that {x, 0) -< (y, 0) if and only if x -< 
y. Without loss of generality, we can assume that the 
number of nonzero Schmidt coefficients of \ip) is n. In 
other words, all Schmidt coefficients of the source state 
are positive. 

For technical simplicity, we apply the above discussions not 
only to normalized but also to unnormalized states. Sometimes 
we shall use a clause such as l |x) can save some entanglement 
lost for the transformation of \ip) to \ip)', and we shall say 
that \ip) can transfer some entanglement into the state |x) 
whenever |x) can do partial entanglement recovery for the 
transformation from \ip) to some unspecified target state \<p) 
to mean that '|x) can do partial entanglement recovery for the 
transformation of \ip) to \<p)'. 

III. Partial entanglement recovery for a 

TRANSFORMATION BETWEEN STATES WITH STRICT 
MAJORIZATION 

In this section, we focus on whether a given auxiliary state 
|x) can do partial entanglement recovery for a transformation 



of \tp) to | ip) such that \rp) is in S°(\ip)). A necessary and 
sufficient condition for such a recovery is presented. Thus, 
a complete characterization of such auxiliary states |x) is 
obtained. 

First, we define the distance between \ip) and \ip) to be 
the Euclidean distance between two n-dimensional ordered 
probability vectors tp* and ip\ i.e, 



\\W-\<p)\\ 



\ 



1 \2 



Before presenting the main result of this section, we prove 
a useful theorem. Assuming that \ip) is in S°(\cp}), we shall 
prove that if £ u (|x)) > 9u(\f)) then a suitable collective 
operation can transform the joint state \ip) ® |x) into another 
joint state \ip) ® \oj) such that \oj) is not 'far from' |x). 
Surprisingly, this result does not depend on which source state 
j^} we have chosen at the beginning. 

Theorem 3.1: Let and \<p) be two states with S 
S°(\ip)). If |x) is an auxiliary state such that l u (\x)) > 
g u (\(p)), then there exists a positive number 5 such that for 
any state \u>) satisfying |||cj) — |x)|| < $, it holds that 

|V) ® |x) ^ \<P) ® 1"). (5) 
Proof. Since l u {\ X }) > gM) and \i>) S S°(\ip)), it 
follows from Lemma 12.31 that 



(6) 



-0 <8> X < f ® X- 



Notice that a small enough perturbation on the right hand 
side of Eq. © will not change the relation '<' since every 
inequality in Eq. (HJ is strict. Thus it is possible to take a 
sufficiently small positive number 6 such that for any state 
\u) satisfying |||o;) — |x)|| < S, the relation tp ® x ~4 f ® w 
holds, which confirms the validity of Eq. (0. With that we 
complete the proof of Theorem 13.11 □ 



The following simple corollary of Theorem 13.11 establishes 
a connection between uniformity indices and partial entangle- 
ment recovery. 

Corollary 3.1: If g u (\ip)) < l u (\x)) < 1> then |x) can do 
partial entanglement recovery for any transformation of \ip) to 
\<p) with |0) G S°(\<p)). 

Intuitively, if the minimal local uniformity of a partially 
entangled pure state |x) is larger than the global uniformity 
of \<p), then the transformation of \ip) to \cp) such that tp <J ip 
can always increase the entanglement degree of |x). 

Example 3.1: Let \tp) and \<p) be two 2x2 states with 
tp^ = (a, 1— a) and <p* = (b, 1 — 6), where i < a < b < 1. The 
goal here is to find a 2x 2 state that can do partial entanglement 
recovery for the transformation of \ip) to \<p). 

Take an auxiliary state \x(p)} with x (p) = {pA ~ P)- It 
is easy to check that \ip) is in S°(\ip)). By Corollary 13.11 if 
|x(p)) satisfies 



9u{\<f)) < lu(\x(p))) < 1, 



(7) 



then |x(p)) can be used to do partial entanglement recovery 
for the transformation \ip) — > \ip). It is easy to see that Eq. (j7]) 
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is equivalent to 



1 



b 1 
< — 



P 



< 1, 



b p 

or i < p < b. The desired state \ui) such that both 
1^) ® |x) — * 1 9?) ® |w) and |w) — > |x) hold can be taken as 
l^- 1 ) — |x(P — e )) w i tn a suitably small positive number e. It 
is obvious that |x) -*» |w) whenever e is positive but small 
enough. □ 

In Example 13. II the condition of p < b means that \x(p)) 
is more entangled than \<p). A simple analysis shows that 
this condition is also necessary to guarantee that \x(p)) does 
partial entanglement recovery for the transformation of \tp) 
to | <p). So we rediscover the main result in [10]: for 2x2- 
dimensional states, the auxiliary state |x) can do nontrivial 
partial entanglement recovery for a transformation with target 
state \p) if and only if \x) is more entangled than \(p). 

Example 3.2: This is a generalization of Example 13. II Let 
\tp) and \p) be two states such that \ip) is in S°(\ip)). Our aim 
here is to decide whether there exists some 2x2 state that 
can do partial entanglement recovery for the transformation of 

W) to \<p)- 

Take an auxiliary state \x(p)) with X (p) = {Pi 1 ~~ P)- By 
Corollary if 



9*(\<P)) < U\X(P))) < 1 



(8) 



then \x(p)) can do partial entanglement recovery for the 
transformation of \tp) to \ip). Moreover, Eq. ([8]) is equivalent 
to 

(9) 



\ <P< 1 



1 



9u{W)Y 

Therefore, the entanglement lost in the transformation of 
\tp) to \tp) can always be partially recovered by a 2 x 2 state 
\x(p)) satisfying Eq. (O. Again, the desired state \u>) such 
that both \ip) (g> |x) — > \ip) (g> \u>) and \u>) — > |%) can be taken 
as \u>) = \x(p— e )) with a suitably small positive number e.D 



In Example 13.21 we show that the entanglement lost in a 
transformation of \ip) to \<p) such that \ip) E S°(\p)) can 
always be partially recovered by a 2 x 2 state \x), and the 
explicit construction of such a state \x) is also presented. This 
is a considerable refinement of Theorem 1 in [12]. We also 
point out that in the proof of Theorem 1 in [12], an important 
constraint on p, i.e., p/3 n < or l u (\x(p))) > g u {\p)), 

is missing, therefore the case (ii) in the proof of Theorem 1 
in [12] is possible if x = s = n and y = t = 0, which makes 
the proof there invalid. 

Corollarv l3.1l onlv provides us with a sufficient condition for 
which \x) can be used to receive some entanglement lost in a 
transformation of \ip) to \<p) with \ip) € S°(\<p)). However, 
this condition is too strong to be satisfied in many cases. 
Nevertheless, the following theorem gives a weaker condition, 
and indeed it provides a complete characterization of states 
|x) that can be used to do partial entanglement recovery for 
a transformation with target state \ip) and source state \ip) in 

Theorem 3.2: Let and \<p) be two n x n states such 
that \if>) is in 5°(|<^)), and let |x) be a k x k auxiliary 



state. Then |x) can do partial entanglement recovery for the 
transformation of to | ip) if and only if one of the following 
three cases holds: 

(i) L u (\x)) — and na > n'(a + 1). Here a and n' are 
the numbers of nonzero Schmidt coefficients of |x) and \ip), 
respectively; 

(ii) g u {\cp)) < L u {\ x )) < 1; 

(hi) L„(|x» = 9u(\v)) and ^ - (x' Bm /C)K Here X [ is 
a segment of x^ with only two distinct components, C is a 
normalization factor, and m > 1. 

Moreover, if none of the above cases holds, then |x) cannot 
do partial entanglement recovery for any transformation of | ip) 
to \ip) such that is in S(\(p)). 

Proof. See Appendix B. □ 

To better understand the above theorem, we give the fol- 
lowing remarks: 

1) The case that L u (\x)) = 1 is not included in cases (i), 
(ii), and (iii). Hence a simple corollary of Theorem 13.21 
is that a maximally entangled state cannot be used to do 
partial entanglement recovery. This is reasonable since 
for such a state |x), there does not exist a k x k state 
\ui) which is more entangled than |x). 

2) The case that L u (\x)) = is slightly different from 
the above case and is more interesting. Although |x) 
is a maximally entangled state in a state space of lower 
dimension ax a with a < k, it is only partially entangled 
when it is considered as a k x k state. Hence it is 
still possible to transform ® |x) into another state 
\tp) eg) |u>), where \u>) is more entangled than |x). Case 
(i) shows that the necessary and sufficient condition is 
that the Schmidt numbers of |x) and \<p) satisfy a simple 
inequality na > n'(a + 1). In some sense, the solution 
in this case explains why the dimensions of the states 
need to be fixed. 

3) Case (ii) means that if |x) is partially entangled and 
the maximal local uniformity of \x) is larger than the 
global uniformity of \<p), then |x) can be used to save 
some entanglement lost in the transformation of to 
| cp). This case provides a feasible sufficient condition 
for partial entanglement recovery. 

4) Case (iii) is of special interest. It supplies the solution 
at the critical point L u (\x)) = 9u{\f))- As we will see, 
the proof of this case is very complicated. We include 
this case for the following two reasons. First, from the 
aspect of the completeness of the solution. Including 
such a special case enables us to completely solve the 
feasibility of partial entanglement recovery for all ip 
and ip with ip < p. Second, from the special form that 
\p>) should satisfy. A careful observation shows that \<p) 
has only two different Schmidt coefficients and should 
be constructed by repeating a segment of x^ finitely 
many times. In our opinion this provides new insight 
into the process of partial entanglement recovery. In 
addition, in the proof of this case we have extensively 
employed the techniques introduced in the present paper 
and the properties of majorization. Hopefully, these 
proof techniques will be useful in solving other problems 



in quantum entanglement theory. 

In sum, Theorem 13.21 provides a necessary and sufficient 
condition under which \\) can do partial entanglement recov- 
ery for some transformation with the target state | ip) . Therefore 
it can be treated as a basic result about partial entanglement 
recovery. In view of Theorem 13.21 it seems reasonable to use 
maximal local uniformity to describe the partial entanglement 
recovery power of an auxiliary state. 

It is worth noting that in the above theorem, the choice of 
has some free degree. That is, if \x) can be used to do 
partial entanglement recovery for a transformation of to 
\(f) such that \ip) is in S°(\tp)), then for any \ip') S S°(\(p)), 
\x) can also recover entanglement lost in the transformation 
of W) to \<p). 

Theorem 13.21 has many interesting corollaries. We only 
consider the following one where the auxiliary state \x) is 
2 x 2-dimensional. 

Corollary 3.2: If |%) and \ip) are two partially entangled 

states with x — (P> 1 ~ P) an d f' 1 — {Pit ■ ■ • > Pn)< men 
\x) can be used to do partial entanglement recovery for the 
transformation of to \ip) such that is in S°(\ip)) if and 
only if one of the following three cases holds: 

(i) p = 1 and n > 2n', where n' is the number of nonzero 
components of ip; 

(iii) p = and \ip) has a special form such that ip- 1 = 

((p, 1 -p)® fc /fc) J - for some k > 1. 

Moreover, if none of (i)-(iii) is satisfied, then \x) cannot do 
partial entanglement recovery for any transformation of to 
\ip) such that \tp) is in S(\ip)). 

The most interesting part of the Corollary 13.21 is case (iii). 
The following example demonstrates this point. 

Example 3.3: Let \x), \<p'), W), and \ip"') be four states 
with x^ = (p, 1 — p)i f = {p,P, 1 — Pi 1 — p)/2, f''^ = 
(p,p, 1 - p)/(l + p), and <p"'± = (p, 1 - p, 1 - p)/(2 - p), 
where -| < p < 1. Obviously, 

Lu(\x)) = 9u{W')) = 9u(\<p")) = Su(|^"» = — • 

By Corollary 13.21 it is easy to see that \x) can do partial 
entanglement recovery for any transformation of \ijj) to \ip') 
with \ip) e S°(\(p')) since (p /l = ((p, 1 -p)® 2 /2) 1 . 

However, again by the above corollary, \x) cannot recover 
anything for any transformations with the target states \ip") 
or \<p"'). ' " □ 

Until now we only deal with the transformations of to 
\ip) such that tp is strictly majorized by cp. What about the 
other cases? In next section, we will prove two more general 
theorems about partial entanglement recovery where \ip) and 
\(p) only need to satisfy the non-strict majorization relation 
ip -< (p. 

IV. Partial entanglement recovery for a general 

TRANSFORMATION 

In this section we deal with partial entanglement recovery 
for a class of more general transformations. Before proceeding 



to the main results, it will be helpful to introduce some 
notations. Let x and y be two finite dimensional vectors. 
We write x C y or y □ x if x[ < y[ and x\ > y\. 
Roughly speaking, x C y means that the values of the extreme 
components of x are strictly bounded by those of y. We use 
the formal expression ^7 □ ^7 as a convenient rewriting of 
x' (g) y" □ x" ® y'. 

For simplicity, in this section we only deal with vectors 
that are already in non-increasing order. That is, for a finite 
dimensional vector x, we assume that x = x^-. 

We now introduce the following concept. 

Definition 4.1: A decomposition of a vector x is a sequence 
of vectors x , . . . , x m satisfying 

(i) each of these vectors has dimension at least one, i.e., 
dim(x I ) > 1; and 

(ii) x is the direct sum of these vectors, i.e., x = 
(x 1 , . . . , x m ), or simply, x = ©™ 1 x m . 

The following simple lemma provides a special decompo- 
sition of two vectors x and y such that x -< y. 

Lemma 4.1: If x and y are vectors satisfying x -< y, then x 
and y can be uniquely decomposed as x — (x , . . . , x m ) and 
y = (y 1 , . . . , y m ) such that 

(i) x l < y % or x l — y % for each i = 1, . . . , m; and 

(ii) there does not exist an index i such that x l — y % and 
x l+1 = y l+1 hold simultaneously. 

Proof. The proof is simple, and the details are omitted. □ 

The decompositions of x and y in Lemma 14.11 are called 
the normal decompositions of x and y. 

Motivated by Lemma 14.11 we shall define two index sets 
I x<y and D XtV for any vectors x and y satisfying x -< y. 
Suppose that x and y have normal decompositions as in 
Lemma [4T1 Then we define 

Ix.y = {i '■ x l — y l and 1 < i < m} 

and 

D x ,y — {i : x l <y l and 1 < i < m}. 
It is obvious that 

Ix, y n L's.y = and I x<y U ZJ^.y = {1, . . . ,m}. 

One can easily check that x < y is equivalent to 7^ = 
and fl^j = {1}. 

In what follows, we only consider the auxiliary state with 
positive Schmidt coefficients, as our major purpose here is to 
find the states that can do partial entanglement recovery for a 
given transformation. For simplicity, the maximally entangled 
state is also not considered. 

With these preliminaries, we present one of the main results 
in this section, which gives a sufficient condition under which 
\x) can do partial entanglement recovery for a transformation 
of \tP) to \<p). 

Theorem 4.1: Let and \ip) be two states with normal 
decompositions t/j = (ip 1 , . . . ,ip m ) and ip — (ip 1 , . . . , ip m ) 
such that -0 is majorized by ip , and let \x) be an auxiliary 
state with a similar decomposition to and \ip), say, x = 

U 1 \'";. if 

X % (p l 

— 7 □ — r, for all i G Iii, v and j € D^i, v (10) 

X^ (p3 
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and 

min{Z u (|x 1 )) : 1 < i < m} > ma-x{g u (\tp t )) : i £ A/>,¥>}> 

(11) 

then \\) can do partial entanglement recovery for the trans- 
formation of to \<p). 

Moreover, if |%) satisfies Eqs. (TTOb and ( fTTt . then there 
exists a positive number S such that for any state \u>) with 
a decomposition w = (w 1 , . . . , w m ) satisfying 

(i) w' = x l f° r an Y * € J^.V' an d 



(ii) ||a/ — x l \\ < 8 and J2 ^ = £ X* f° r any 



the 



transformation of |^) ® |x) to |<p) ® |w) can be realized with 
certainty by LOCC, i.e., \<tp) ® |x) -> |y) ® 

Proof. See Appendix C. □ 

Here we give some remarks: 

1) In the above theorem whether \x) can save entanglement 
lost for the transformation of \ip) to \tp) does not directly 
depend on the choice of the source state It only 
depends on the decomposition of the target state \tp) and 
the index sets I^ >(p and D^ jlp . For this reason, in what 
follows, it is not necessary to specify the source state 
|?/>) clearly. We only need to give a decomposition of tp 
and two index sets / and D. Based on these conditions, 
we can identify a class of auxiliary states \x) that can 
do partial entanglement recovery for any transformation 
of \ip) to \tp) with source state \ip) satisfying J^ tlp — I 
and = D. 

2) When \x) can be used to do partial entanglement re- 
covery for the transformation of \ip) to \tp), the more 
entangled state \u>) generated from \x) is also explicitly 
given by the above theorem. It should be noted that in 
general the resulting state \w) is determined by the states 
|x), IV*) and \tp) together, although the choice of \x) 
doesn't depend on the source state \tjj). In other words, 
sometimes there may not exist a universal state \oj) in 
the sense that \ip') ® |x) — > \p) ® \oj), \u) — » \x) and 
|x) \w) hold for all states \ip') with D^' tV — D^, tip 
and Iij)\ip = Iiii,ip- 

We now examine some special cases of Theorem 14.11 The 
first special case is that both D and I are singletons. 

Corollary 4.1: Let \tp) and |x) be two states with decom- 
positions tp = (p l ,tp 2 ) and x = (xSx 2 ), and let / = {1} 
and D = {2}. If \ip) and |x) satisfy 



X 



□ 
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and 



i{lu(\x 1 )),U\x 2 ))}>9u(\<p 2 )), 



(12) 



(13) 



then |x) can do partial entanglement recovery for any transfor- 
mation of to \<p) with source state \t/j) such that I^ tV = I 
and D^ v = D. 

A corresponding result for the dual of case / = {2} and 
D = {1} can be obtained by exchanging x 1 with x 2 and cp 1 
with ip 2 in Eq. (fT2l and Eq. (13[ , respectively. 

For the sake of convenience, for two n-dimensional vectors 
x and y with x -< y, we define A X)J/ as the set of all indices 



m such that the inequality in Eq. ([TJ holds with an equality, 
i.e., 



A 



x,y 



{m : e m (x) = e m (y) and 1 < m < n — 1}. 



Note that 1 G A Xiy is equivalent to x\ = y% and n — 1 e A^, jV3 
is equivalent to x n — y n . 

Now we present two examples to illustrate the use of 
Corollary [4~T1 

Example 4.1: Let \ip) and \ip) be two nxn states such that 
\tp) is in S(\p)) (n > 2). Assume that A^ ^, = {1}. We hope 
to find an auxiliary state |x) with the minimal dimension to 
do partial entanglement recovery for the transformation of \tp) 
to \p). 

To be more specific, let tp = ((3±, . . . ,(3 n ). Since A^, jV = 
{1}, it is obvious that tp has a normal decomposition <p = 
(tp 1 , tp 2 ), where tp 1 = (Pi) and p 2 — (fi 2 , ■ • ■ , Pn)- Moreover, 
Iil>,u> = {1} and D^ tp — {2}. Take an auxiliary state |x) 
with x = (x\x 2 ), where x 1 = (71), X 2 = (72,73), and 
7i > 72 > 73 > 0. 

By Corollary gj] if |x) satisfies Eqs. P-dUl then |x) 
can do partial entanglement recovery for the transformation of 
to \tp). So we have 



71/72 > P1/P2, 

7l/73 < Pl/Pn, 



(14) 
(15) 



and 



73/72 > Pn/fo. (16) 

By Eqs. ( fT4b and ( TToT ), we can take positive numbers A and 



/i such that 



and 



71 =72(l + A)/3i//3 2 



73 = 72(1 + fl)P n /02- 



(17) 



(18) 



Substituting Eqs. ([T7]> and (|T8j into Eq. dT3J yields < A < 
fi. Moreover, the constraint 72 > 73 and Eq. ( fT8l yield /i < 
(02 — Pn)/Pn- 72 is used to make the following normalization 
condition satisfied: 



3 

E7» 



1. 



(19) 



Notice that P2 > P n . One can easily check that such 
(71,72,73) satisfying Eqs. dTTb— (fl~9b is a solution of the sys- 
tem of inequalities defined by Eqs. dT4b— (fT6b. The parameters 
A and fi satisfy < A < /x < (p 2 — Pn)/Pn- 

Thus by Corollary 14.11 the 3x3 auxiliary state |x) can 
do partial entanglement recovery for the transformation 
of \ip) to \tp). Moreover, the state \aj) such that both 
® |x) ~ * \f) ® l w ) and \lu) — > |x) hold can be chosen 
as ui = (71,72 — e,73 + e) with a sufficiently small positive 
number e. □ 

We point out that the existence of such an auxiliary state 
|x) with x = (71,72,73) has been proven in Theorem 3 in 
[12], where |x) is of the form x(p, <?) = (p, 9, 1 ~ P — ?)> 
P><7>1— p — q > 0. However, an important constraint 
on the p and q or the region R, i.e., (1 — p — g)/?2 > ?/3n, is 
missing in [12], Thus an additional case which is not included 
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in Case (i) and Case (ii) in [12] is possible, which makes that 
proof invalid. 

A dual case of Example 14. II is as follows: 

Example 4.2: Let \ip) and \tp) be two nxn states such that 
\ip) is in S(\tp)) (n > 2). Assume that A^, v = {n — 1}. We 
hope to find an auxiliary state |x) with the minimal dimension 
to do partial entanglement recovery for the transformation of 

W) to \<p)- 

To be more specific, let tp = (J3t, . . . , (3 n ). Since A^ tV = 
{n — 1}, it is easy to check that tp has a normal decomposition 
tp = (tp 1 ^ 2 ), where tp 1 = {fl x , . . . ,0 n -i) and tp 2 = (/?„). 
Similarly, let an auxiliary sate |x) have a decomposition \ = 

(x\x 2 ), where X 1 = (71,72), X 2 = (73), and 71 > 72 > 
73 > 0. By Eqs. (fT2l and ( fT3] l again, noticing that 1^ „ = {2} 
and .D^, v = {1}, we have the following system of inequalities: 



and 



73/71 > Pn/Pl, 
73/72 < Pn/Pn-l 



(20) 
(21) 



and 



72/71 > Pn-l/Pl- (22) 
By using a similar argument as in Example 14.11 we can take 

72 =7i(l + M)/3„-i//?i (23) 

and 

73 = 7i(l + A)/V/3i, (24) 

where < A < p < — /3„_i)//3„_i. (A < /i is deduced 
by substituting Eqs. d33]l and (fJHi into Eq. (EB, M < - 
f3 n -i )//3 n _i comes from Eq. (l23l and 71 > 72). 71 is taken 
to validate the following normalization condition 



3 

i=i 



1. 



(25) 



Since /3i > one can easily check that such a state 

|x) with x = (71,72,73) is a solution of the inequalities 
system defined by Eqs. (|20]>-(|221. Thus by Corollary 14. 1 1 |x) 
can do partial entanglement recovery for the transformation 
of \ip) to \tp). Again, the desired state \u>) such that both 
® \x) ~ ¥ \<p) ® l w ) and l w ) — * Ix) hold can be chosen as 
u> = (71 — e, 72 + e, 73) with a suitably small positive number 
e. □ 

If one of the cases A^ i¥ , = {1} or „ = {n— 1} occurs, 
we can always use a 3 x 3 state \x) to partially recover 
entanglement lost in the transformation of to \<p). The 
explicit construction of such \\) has also been presented in 
the above examples. 

The following corollary is another important special case of 
Theorem 14.11 

Corollary 4.2: Let \<p) and \x) be two states with ip — 
(tp\<p 2 ,<p 3 ) and X - (x\x 2 ,X 3 )- / - {1,3} and D = {2}. 
If |x) and \<p) satisfy 



X 2 v 21 

A_ □ __. 

X 2 V 2 ' 



(26) 



(27) 



min{^(|x 1 )),Zu(|x 2 )),^(lx 3 ))}>5«(^ 2 )), 



(28) 



then |x) can do partial entanglement recovery for any trans- 
formation of I tp) to \ip) with p = I and = D. 

A very interesting application of the above corollary is the 
following: 

Example 4.3: Let and \<p) be two nxn states such 
that \ip) is in S(\ip)) (n > 3). Assume A^ >ip = {l,n — 1}. 
Our purpose here is to find an auxiliary state to do partial 
entanglement recovery for the transformation of \ip) to \ip). 

Take a 4x 4 auxiliary state |x) with x = (71, 72, 73, 74)- Let 
us decompose tp and x, respectively, into tp — (tp 1 , tp 2 , tp 3 ) and 
X = (x\x 2 ,X 3 ), where tp 1 = (fix), ip 2 = (fo, . . . , /3„_i), 
f 3 = (/?n). X 1 = (7i), X 2 = (72,73), X 3 = (74), and 
7i > 72 > 73 > 74 > 0. Since A^ <(p = {l,n - 1}, it is 
easy to check that I$,u> — {1,3} and D^ iV — {2}. Thus by 
Corollary 14.21 |x) can do partial entanglement recovery for the 
transformation of \ip) to \tp) if Eqs. (|26]>-(|28]> hold. A routine 
calculation leads to the following solution of Eqs. (|26T> — (f28l>: 

71 =72(l + A)A/#2, 
73 = 72(1 + ??)(! +ll)Pn-l/fa, 

74 = 72(1 + A 1 )/V/32, 

where 1 + A < (1 + rj)(l + /1) < /3 2 //3„-i, A, p, and r; 
are all positive real numbers, and 72 is used to validate the 
normalization condition 
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E 

i=l 



7i 



1. 



So such an auxiliary state |x) for partial entanglement recovery 
always exists. 

Ultimately, to partially recover entanglement lost in the 
transformation of \ip) to \tp), it is sufficient to use an auxiliary 
state |x) with dimension 4x4. Again, the more entangled 
state \uj) generated from |x) after the recovery process can be 
chosen as uj = (71 , 72 — e, 73 + e, 74), where e is a sufficiently 
small positive number. □ 

In [12], it is proven that any 3x3 state cannot be used 
to partially recover entanglement lost in the transformation of 
to \tp) with ip -< tp and A^. v = {1, n — 1}. By the above 
example, we are able to show that 4x4 auxiliary states are 
necessary and sufficient to do partial entanglement recovery 
for this special case. 

In practice, we hope that the dimension of the auxiliary 
state |x) is as small as possible. In Theorem 14.11 if there are 
two successive integers i and i + 1 both contained in D^ !V (in 
Iif, t(p this case cannot happen), we in fact can combine x 1 with 
X 1+1 to reduce the dimension of \. So a careful investigation 
of the structure of is necessary. 

Let us see a simple example. Suppose that for states \ip) and 
\tp), = {1, 4, 7, 12} and = {2, 3, 5, 6, 8, 9, 10, 11}. 

By the construction in Theorem 14.11 we should use an 
auxiliary state |x) with \ = (x , ■ ■ ■ , X 12 ), where each 
X* (i € D^ tip ) has dimension at least 2. Thus the vector 
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20. If we and 



X has dimension at least | 

I + 2 |A/,, V | 

combine the successive integers in D^ tV together, we have 
D'^ = {{2, 3}, {5, 6}, {8, 9, 10, 11}}, and the dimension of 

10. 

) and \tp) are two states 
(ip 1 , . . . , ip m ) and tp = 



2 id:, 



X is reduced to \I^ >lp \ - r ~^ iV> \ 
More formally, suppose that 
with normal decompositions ip 
(ip 1 ,...,?™). Let 



where 



={h,... ,k p }, 1 < p < 



= fco < fci < • • • < k p < k p+ i = m + 1. 



We define 



where 



D' = {Di : A + and < % < p}, 



Di = {s:h + l<s< k l+1 - 1}. 

The constraint Di ^ in the definition of D'^ needs a careful 
explanation. For any 1 < i < p, we have that ki £ 1^ „ implies 
ki + 1 ^ Ii/>,<p- So Di 7^ in this case. However, if kx = 1 
or k p = m then we have Do = or D p = 0, respectively. 
To avoid these two trivial cases, the constraint Di ^ is 
necessary. In particular, if I^ >v — then D'^ — {D^ i(p } = 
{{1, . . . , m}}. For the sake of convenience, we also define 

4>,v = {{*} : 1 e h,<f}- 

In the following discussions, we shall use the elements of 
1^ and DL as indices. We define the natural order of the 
elements in I^ >(p U D^ as 

D < {fci} < Di < {fc 2 } < ■ ■ ■ < {k p } < D p , 

where we assume that any term which doesn't exist should 
be omitted automatically without affecting the orders of other 
terms. 

Suppose that J is a finite set of integers. We use the 
notations max J and min J to denote the maximal and the 
minimal elements of J, respectively. For any real function /(.) 
defined on J, the expression arg min fc6J /(fc) denotes the 
index i E J such that f(i) — minfegj/(fc) (here we assume 
that there is a unique i of J that can attain the minimum). 

Now we can present another condition for the existence of 
partial entanglement recovery, which complements Theorem 
EJJ 

Theorem 4.2: Let and \ip) be two states with normal de- 
compositions ip = (ip 1 , . . . ,ip m ) and ip = (ip 1 , . . . , (p m ) such 
that %jj is majorized by tp, and let \x) be an auxiliary state with 



a decomposition x = (x D ":X^ kl \x I 



.X 



htl', UD' X 



if 



(29) 



for all i e J € D'^ and Jj = arg min fce 7 |i - k\, 



min{/„(|x J )) : J G DL } > max{ 5u (|^)) : i e D^ v }, 



min{I„(|xW» : i € W > max \J {g u {\^ J )), 

(31) 

then \x) can do partial entanglement recovery for the trans- 
formation of \t/j) to \<p). 

Moreover, if \x) satisfies Eqs. (|29b-(f3Tb, then there exists 
a positive number 5 such that for any state \u>) with a 
decomposition uo = © ie /' u l satisfying 



(i) u> 1 — x l for an y i ^ I'll 



and 



(ii) || uj 1 - xil < S and E to 1 = X % for any i e D'^, the 
transformation of \ip) ® |x) to |</?) ® |a->) can be realized with 
certainty by LOCC, i.e., \%p) ® \x) -* \<p) ® 
Proof. The proof is similar to Theorem 14.11 and we omit the 
details. □ 

The key idea in the above theorem is to let all vectors 
between (p ki+1 and (p ki + 1 ~ 1 correspond to a single x Di ■ This 
reduces the dimension of \ efficiently. 

An interesting special case of Theorem 14.21 is when the 
majorization ip -< ip splits into m strict majorizations: ip' 1 <jLp l . 
We state this result in the following: 

Corollary 4.3: Let \ip) and \<p) be two states such that \ip) is 
in S(\(p)). Suppose that ip and tp have normal decompositions 
ip = (ip 1 , ip m ) andip^(p\..., tp" 1 ), and let = 0. 
If |x) is an auxiliary state such that 

lu(\x)) > max{g u (\tp i )) : 1 < % < m}, 

then \x) can do partial entanglement recovery for the trans- 
formation of |^) to \tp). 

Proof. In fact, in this special case, D^ = {D^, i¥ ,} = 
{{1, . . . , m}}, IL = 0. Thus, by Theorem 14.21 to do partial 
entanglement recovery, the only non-trivial condition that \x) 
should satisfy is Eq. (f30b . which is exactly the assumption of 
the present corollary. □ 

It is easy to check that in the above corollary \x) can be 
chosen as a 2 x 2 state. However, by Theorem 14.11 we can 
only find a state \x) of dimension at least 2m x 2m. 

By summarizing Theorems 13.21 14.11 and 14.21 we have the 
following: 

Theorem 4.3: Suppose that \ip) and \tp) are two nxn states 
such that ip -i, tp. We can always find an auxiliary state \x) to 
do partial entanglement recovery for the transformation of \ip) 
to \ip), where the dimension of x is between 2x2 and nxn. 
Moreover, such a state |x) can only depend on the target state 
\tp) and the presence of equalities in the majorization ip -< tp. 

Proof. The proof follows immediately from Theorems 13.21 
RTTl and |431 □ 



The upper bound nxn cannot always be reduced to (n — 
1) x (n— 1). We have seen in Example 03] that when n = 4, an 
auxiliary state \x) of dimension 4 x 4 is needed to do partial 
entanglement recovery for the transformation of \ip) to \tp) 



(30) such that A 



{1,3}. 
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We conclude this section by giving an example to illustrate 
the use of Theorem 14.21 This example is taken from [13]. 

Example 4.4: Let and \ip) be two n x n states such 
that \tp) is in S(\cp)) {n > 6). Assume that A^ i¥> = {2,3,5}. 
The goal here is to find an auxiliary state |x) to do partial 
entanglement recovery for the transformation of \tj}) to \<p). 

To be specific, let ip = (fa, . . . ,fa). It is easy to check 
that cp has a normal decomposition <p — (ip 1 , p 2 , p 3 , (p 4 ) , 
where <p l = (0i,0 2 ), ^ 2 = (0 3 ), V 3 = (fa, fa), and 
p 4 = (fa, . . . ,0„). Also, 7^ = {2} and D^ v = {1,3,4}. 
So = {{2}} and 0^ - {{1}, {3, 4}}. 

Take a 5 x 5 auxiliary state |x) with x — 

(X {1 \X {2 \X {ZA} \ where X {1} = (71,7a), X {2} - (7s), 
x {3,4} = ( 74j75 ) ? and 7i > • • • > 75 > 0. By Theorem l4~2l 
Eq. d29j yields 



and 



x {2} _ ^ 

x {3,4} J ^3 ■ 



(32) 



(33) 



Eq. d30]l yields 

mm{l u (\ X W)),l u (\ X ^))} > max{ 5u (|^)), 

ff„(|^)),ff u (|^ 4 ))}. 

(34) 

Eq. (1311 yields 

Z U (|X {2} » > max{ 5u (|^)), 3u (|^)), 5u (|^))}, (35) 

which is automatically satisfied since 1 U {\X {2} )) = 1 while the 
right hand side of Eq. ( T35l l is strictly less than 1. 
More explicitly, we have 



and 



and 



73 ^ 03 , 73 ^ 03 
71 Pi 72 P2 



73 03 - 73 ^ 03 

— > — and — < — . 

74 04 75 05 



. r 72 75s ^ r 02 05 071-, 

mm{ — , — > max{ — , — , — }. 

71 74 01 04 06 



(36) 



(37) 



(38) 



With a routine calculation one can check that 

0i 

7i = M7375-, 

03 

72 = (l + »?)73-A 

P3 

, 04 

74 = h^ — , 

fa 



75 = (1 + A)73 



is a solution of the system of inequalities defined by Eqs. 
(f36T> — (f38l>. where 73 > is used to satisfy the normalization 
condition: 



The parameters r\, A, /i, and h satisfy 

01-02 _ ^ . 04 " 05 



02 



05 



and 



i=l 



7 = 1- 



/-it \02 ^ ^ • r / -j 1 \0204 /, . X0206 -, 

(1 + V)-r < V < mm{(l + i)-—^, (1 + V)-——, 1}, 

01 0105 010n 

1 ^ 1, - • m 1 \\0105 /! , N0506 

(1 + A)— </i<mm{(l + A)— — ,(l + ry)— — ,1}. 

P4 P2P4 P4Pn 

Notice that (3\ > fa, 04 > 05, an d fa > fa, and such a state 
|x) with x — (71, ■ ■ • 175) always exists. So we have actually 
constructed a class of states |x) with dimension 5x5 that 
can do partial entanglement recovery for the transformation 

of \ip) to \<p). □ 



V. A POLYNOMIAL TIME ALGORITHM FOR PARTIAL 
ENTANGLEMENT RECOVERY 

In this section we study partial entanglement recovery from 
the algorithmic viewpoint. We present a polynomial algorithm 
of time complexity 0(n 2 k 4 ) to decide whether |x) can be 
used to recover some entanglement lost in the transformation 
of |^) to \<p), where n and k are the dimensions of ip and x, 
respectively. 

The key part of Problem 1 is to solve the majorization 
relation V^X -< (f®u>. As argued before, the main difficulty 
here is how to deal with the order of the tensor product ip ® u 
when (j varies. We will develop some techniques to overcome 
this difficulty. Notice that the map from u> to ip eg) uj is an 
affine one. To make our discussions more general and more 
readable, in what follows we consider affine maps instead of 
tensor products. 

To be concise, some concepts are introduced first. 

Definition 5.1: Let / and g be real functions defined on 
TZ m , and let S C 1Z™ . f and g are said to be comparable on 
S if 

(i) Vz e S, f(x) > g(x); or 

(ii) Vx g S, f(x) < g(x). 

Let F be a map from K m to TZ n . We write F(x) = 
(fi(x), ■ ■ ■ , f n (x)), where each /, is a real function defined 
on W n . 

Definition 5.2: F is said to have a fixed order on S if for 
any 1 <i < j < n, fi and fj are comparable on S. 

Suppose F has a fixed order on S, and assume that whether 
fi and fj are comparable on S can be determined in 0(1) 
time. Then there exists a common algorithm which can sort 
the entries of F(x) into non-increasing order for any x G S in 
0(n log 2 n) time. This fact is extremely useful in the following 
discussions. 

If some entries of F are not comparable on S, then by the 
above definition F does not have a fixed order. An important 
question naturally arises: how many different orders can F 
have on 5? 

Definition 5.3: F is said to have at most M different orders 
on S if there exists a decomposition of S, say, 6>i, • • • ,Sm, 
such that 

(i) S = Si U • • • U 5 A /; and 
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(ii) F has a fixed order on each i = 1, • • • , M. 

Now let F be an affine map, F(x) = Ax + b, where A £ 
TZ nxm and b £ TZ n . At first glance, F may have n! different 
orders on lZ m . However, this is not always true. A somewhat 
surprising fact is that the number of different orders of F can 
be dramatically reduced to 0(n 2m ) when m is a constant. 

Lemma 5.1: F has at most 0(n 2m ) different orders on lZ m . 

Proof. For any 1 < i < j < n, the difference of fi(x) and 
fj{x) crosses zero (from positive to negative, or negative to 
positive) if and only if x crosses the hyperplane determined 
by the equation /, (x) — fj (x) = 0, or more precisely, 

m 

Tij = {(xi, • • • ,x m ) : ^(a is - a js )x s + {bi - bj) = 0}. 

s=l 

It should be noted that there are two cases where Ty does not 
define a legal hyperplane. The first case is Ty = and the 
second one is = lZ m . We will exclude these cases since in 
both of them /, and fj remain comparable whatever x varies. 
Denote 

T = {Tij : ^ and ^ K m , l<i<j<n}. 

The number of hyperplanes in T is less than or equal 
to n(n — l)/2. These hyperplanes divide TZ m into at most 
0((n(n - l)/2) m ) = 0(n 2m ) different parts. F has a fixed 
order on each part. With that we complete the proof. □ 

It is obvious that the above lemma holds for any subset of 

n m . 

Lemma 15.11 indicates that we can decompose lZ m into 
0(n 2m ) parts, T>%,- ■ ■ ,T>m, such that on each part, F has a 
fixed order. In practice it is important to construct these parts 
explicitly. To see how this procedure can be done efficiently, 
let us first examine a special case where m = 1. 

Example 5.1: Let F(x) — {a%x — b\, • • • , a n x — b n ), where 
x > 0. For simplicity assume a, 7^ dj, bi 7^ bj, a,, 6, > for 
any 1 < i < j < n. 

By Lemma I3TT1 F has at most 0(n 2 ) different orders when 
x varies as a non-negative number. In what follows we will 
show how to determine these orders explicitly. 

Step 1. For each 1 < i < j < n, solve equation a,iX — bi = 
cijX—bj. The solution is given by #y = (bi—bj)/ (a — a,j). Let 
r = {Oij : 1 < i < j < n} U {0}. The number of elements of 
r is denoted by M. It is easy to see that M < n(n — l)/2 + l. 

Step 2. Sort the elements in V into non-decreasing order, 
say = c < c\ < ■ ■ ■ < c M -i- 

Step 3. Construct a sequence of intervals: T>\ = [cq, c{\, ■ ■ •, 
T>m = [cm -1, +00). 

It is clear that F has a fixed order on each interval. The 
above procedure is completed in 0(n 2 ) + 0(M\og 2 M) + 
O(M) = 0(n 2 log 2 n) time. 

It is notable that the leftmost interval [cq, c{\ can be located 
in 0(n 2 ) time. This fact will be useful in the following 
discussions. □ 

To deal with the general case, we need a lemma in com- 
putational geometry. Let 7i be a set of n hyperplanes in lZ d 
with d > 1. Then H divides lZ d into 0(n d ) parts with pairwise 



disjoint interiors. We call the set of these parts a c?-arrangement 
of fi. A celebrated result in computational geometry shows 
that the d-arrangement of 7i can be enumerated efficiently 
[24]. 

Lemma 5.2: The <i-arrangement of n hyperplanes may be 
computed in time 0(n d ). 

Employing Lemma 15.21 we can easily see that the above 
decomposition T>i, ■ ■ ■ , T>u can be computed in 0(n 2m ) time 
in the case that m > 1. 

With the aid of Lemma 15.11 we are able to solve a 
majorization inequality of the form c -< Ax + b by using 
linear programming methods. 

Lemma 5.3: The majorization inequality c -< F(x) can be 
solved in 0(n 2m+1 log 2 n) time, where m is treated as a 
constant. 

Proof By Lemma IBTTI F has at most M = 0(n 2m ) different 
orders on lZ m . Let us decompose lZ m into M parts and 
enumerate them as T>i,---,T>m- On each part T>i, F has 
a fixed order. This procedure needs time 0{n 2m ). In what 
follows we will show on each part, the majorization inequality 
c -< F(x) can be solved in 0{n\og 2 n) time by using 
standard methods of linear programming. Hence we obtain an 
algorithm with time complexity 0(n 2m ) + 0(n 2m n log 2 n) — 
0(n 2m+1 log 2 n) to solve the desired majorization inequality 
on K m . 

Let us concentrate on a specific An algorithm to solve 
the majorization inequality on T>i is as follows: 

Step 1. Sort c and F(x) into non-increasing order, re- 
spectively. Assume = (c' 1 !,-",^"') and F^(x) = 
(oW^-.oWfi)). 

Step 2. Transform the majorization inequality c -< F(x) into 
the following linear system of inequalities: 

1 1 

J2c (s) < ^]a (s) (x), l<l<n, (39) 

8=1 8=1 

with equality holding when I = n. 

Step 3. Solve the system of inequalities in Eq. d39l using 
standard techniques of linear programming. 

Now let us calculate the time complexity of each step. It is 
obvious that c can be sorted non-increasingly in 0(nlog 2 n) 
time. Since F has a fixed order on T>i, F(x) can also be 
sorted into non-increasing order in 0(nlog 2 n) time. So Step 
1 can be completed in 0(n log 2 n) time. To figure out the time 
complexity of Step 2, we need the following simple fact: the 
linear transform of [y x , ■ ■ ■ , y n ) to (y 1 , y 1 + y 2 , ■ ■ ■ , yi + y 2 + 
• • • + Un) needs only 0(n) time. So Step 2 needs 0(nm) = 
0(n) time. The time complexity of Step 3 needs a careful 
analysis. By applying the well-known Karmarkar's algorithm 
in the theory of linear programming [25] directly, Step 3 needs 
0(n 3 5 ) time. However, in [26], it has been shown that linear 
programming can be solved in linear time 0(n) when the 
dimension of variable x is fixed. Hence the total time to solve 
c -n F(x) on Di is 

0(n log 2 n) + 0(n) + 0(n) = 0(n log 2 n). 

With that we complete the proof of Lemma 15.31 □ 
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Now we are able to present our algorithms about partial 
entanglement recovery. The first algorithm to solve Problem 1 
is a direct consequence of Lemma 15.31 

Theorem 5.1: Problem 1 is solvable in 0(n 2fe_1 log 2 n) 
time, where k is treated as a constant. 

Proof. The key here is to solve the majorization inequality 
■0 ® X -4 <p ® us. Notice that when ip is fixed, the map 
from u> to ip® is an affine one. So Lemma 15.31 works. A 
subtle point here is that w is a fc-dimensional probability 
vector and has only k — 1 independent parameters. In 
addition, the relations ui -< x an d X^ 7^ w can easily 
be cast into linear constraints of to. The total number 
of these constraints is at most 0(kl) = 0(1) when 
k is a constant. Hence the time complexity is in fact 
0({nk) 2( - k -V +1 log 2 n) + 0(1) = O^ 2 * 8 " 1 log 2 n). □ 

The main advantage of the above algorithm is that it can 
determine all the resulting states \uS) in the process of partial 
entanglement recovery. However, this algorithm is efficient 
only when k is treated as a constant. If k varies freely, it will 
turn into exponential time complexity and cannot be efficient 
anymore. To further reduce the time complexity, some lemmas 
are necessary. 

Let = (7i) • ' ' ) 7*0- For the sake of convenience, we 
assume all k entries of x are distinct. The general case can 
be considered similarly by using the compact form of x- For 
each 1 < i < j < k and e > we introduce the following 
vector: 

x(M» = (7i>--- ) 7i - e !---,7j +<V,7fc)- (40) 

To keep the order of x(hj, e ) fixed when e varies, the 
constraints 

7i " e > 7i+i and 7i-i > Jj + e 

should be satisfied. Let Sij be (7^ — 7, + i)/2 if j = i+ 1, and 
be min{7i — — jj} otherwise. Then e 6 [0, 6ij]. 

The following two lemmas exhibit some interesting 
properties of the solutions of Problem 1. Interestingly, the 
first lemma shows that we only need to consider the solution 
\ui) with the Schmidt coefficient vector of a special form 
given in Eq. ( 140) . 

Lemma 5.4: Problem 1 has a solution if and only if there 
exist 1 < i < j < k and e 6 (0, Sij] such that ip <S> x ~< 

V®x(i,3,z)- 

Proof. Sufficiency: Suppose such i, j, and e do exist. It is 
easy to verify x(hj> £ ) X an d X 1 7^ X (hj> £ )- These facts, 
together with the hypothesis ip (g) x ~< V ® x(hj> e X indicate 
that x(«,j, e) is a solution of Problem 1. 

Necessity: Assume Problem 1 has a solution \oj). Then we 
have ifj®x^ L P® UJ > UJ ^X> an d ^ ^ X" ■ The existence of 
i, j, and e such that tp (g> (p -< ip ® x(z, j, e) follows directly 
from the following two facts: 

(a) If a; -< x an d X^ 7^ w ^ tnen there exist 1 < i < j < k 
and e € (0, Sij] such that w -< x{hh e ) ~< X- This is a direct 
consequence of Lemma ITTTI in Appendix B. 

(b) Any state \x') such that u> -< x' -< x a nd 7^ x'^ is 
also a solution of Problem 1. This follows directly from our 



formulation of Problem 1 . □ 

Lemma 5.5: If \x(hj> e o)) is a solution of Problem 1, then 
for any < e < eq, \x(hj, e )) is a lso a solution. 

Proof. Immediately from the formulation of Problem 1 and 
Eq. <|40|). □ 

We are now in a position to state the main result of this 
section, the promised algorithm of time complexity 0(n 2 k 4 ). 

Theorem 5.2: Problem 1 is solvable in 0(n 2 k A ) time. 

Proof. By Lemma 15.41 we only need to consider the 
following problem: for each specific pair (i,j) such that 
1 < i < j < k, decide whether there exists e € (0, S^] such 
that ip ® x -< <p ® x{i 1 j, e )- In what follows we show that this 
problem can be solved in 0(n 2 k 2 ) time. Then by enumerating 
all possible pairs of (i,j), we get an 0(k(k — l)/2n 2 k 2 ) = 
0(n 2 k A ) time algorithm to solve Problem 1. 

Let us begin with two specific indices i and j. By Lemma 
15.11 the number of the different orders of (p ® xihf e ) is a t 
most 0((nk) 2 ) when e varies in [0, Sij]. With Lemma IB31 in 
mind, it is enough to consider one special order among them. 
More precisely, suppose the interval [0, Sij] is divided into M 
parts (intervals), namely, 

V-y = [co,Ci],T>2 = [C1,C 2 ],--- ,T> M = [c M -l,C M ], 

where = cq < c\ < ■ ■ ■ < cm = Sij, and M — 0((nk) 2 ). 
On each interval (p®x(h f e ) nas a fixed order. By lemma [531 
if \x(h ji e o)) is a solution of Problem 1, then any \x(i,j,e)) 
such that < e < min{eo, ci} is also a solution. So we need 
only to consider the leftmost interval T)\ . Our algorithm goes 
as follows: 

Step 1: Find c\\ 

Step 2: Sort ^ ® x and f ® x(h3> e ) i nto non-increasing 
order, respectively, where e € [0, ci]; 

Step 3: Solve the system of inequalities induced by the 
majorization relation ip <g> \ f ® x(*j j; e )- 

Step 4: Output: if a solution of e > is obtained in Step 3, 
then Problem 1 has a solution \x(hj, e ))> otherwise Problem 
1 does not has a solution of the form \x(i,j,e)) for fixed i 
and j, and e G (0, <%]. 

Step 1 requires that we search for the smallest positive 
elements among M items, which requires 0(M) = 0((nk) 2 ) 
time (see also Example l5.ll ). Step 2 needs 0(nk log 2 nk) time. 
Step 3 merely needs 0(nk) time since there is only a single 
parameter e. Step 4 only needs 0(1) time. In sum, only 

0((nk) 2 ) + 0(nk log 2 nk) + 0(nk) + 0(1) = 0((nk) 2 ) 

time is required. □ 

In view of Theorem 15.21 we can say that Problem 1 is 
efficiently solvable. It also suggests that we can study the 
process of partial entanglement recovery using algorithmic 
methods. 

To conclude our discussions about Problem 1, we would 
like to address an important issue for further study. In almost 
all the results we obtained so far, we are only concerned 
with the feasibility of partial entanglement recovery, while the 
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efficiency of this process has not been touched yet. These 
results are of limited use in practice, when we hope to 
minimize entanglement lost in LOCC transformations. In other 
words, we require the resulting state \w) to be not only more 
entangled than |x), but also an "optimal" one that we can 
achieve in this process. Using entropy of entanglement as 
a measure, we suggest the following optimization problem. 
We also note that some aspects of the efficiency of partial 
entanglement recovery have been discussed in [14]. 

Open problem: Given a triple of states (\ip) , \ ip) , \ \)) such 
that ip -< (p, let ft — {uj : ip (g) \ -< <p ® oj and 10 < x}- 
Maximize E(\oj)), subject to u> G Q. 

In the above problem we remove the constraint x 1 
u>*. This makes £1 compact. Thus the continuous function 
E(\u)) can attain its maximum on f2. Suppose \uo) is one 
of the states attaining the maximum. Noticing that the entropy 
of entanglement decreases under majorization, we have the 
following simple relation 

E(\ X )) < EQwo)) < E(\iP)) - E{\p)) + E(\ X )), 

where the first inequality is from uiq -< x, and the second 
inequality is from ip ® \ -< p (gi ojq and the additivity of 
entropy of entanglement. The first inequality is an equality 
if and only if u* — x 1 f° r an Y w G fi, i.e., |x) cannot do 
partial entanglement recovery for the transformation of to 
| cp) . The second inequality is an equality if and only if (ip (g> 
xY = (<p <S) ■ Theorem 15.21 in fact provides a polynomial 
time algorithm to determine whether the first inequality holds 
strictly. How to design efficient algorithms to find the optimal 
state \ujq) seems to be a challenging and worthwhile problem. 

VI. Some applications 

In this section, we establish some interesting connections of 
partial entanglement recovery to the generation of maximally 
entangled states, quantum catalysis, mutual catalysis, and 
multiple-copy entanglement transformation. 

A. How to obtain maximally entangled states by using partial 
entanglement recovery 

Maximally entangled states play a crucial role in many strik- 
ing applications of quantum entanglement such as quantum 
superdense coding [2] and quantum teleportation [3]. It is 
very important to generate such states in practical information 
processing. Under the constraint of LOCC, a natural way 
to obtain a maximally entangled state is to concentrate a 
large number of partially entangled states [5], However, such 
a concentrating protocol involves infinitely many copies of 
the source state while in practice only finitely many copies 
can be available. One can find various deterministic protocols 
based on Nielsen's theorem [6] and probabilistic protocols 
based on Vidal's theorem [17] (see also [16]). It has been 
shown that two 2x2 partially entangled states sometimes can 
be concentrated into an EPR pair deterministically [10]. An 
extensive generalization of such a deterministic concentration 
protocol was presented in [11], where the maximal number 
of Bell states that can be concentrated from a finite number 



of partially entangled states was derived. In what follows, we 
consider deterministic transformations only. 

The following theorem shows that almost all deterministic 
entanglement transformations can concentrate a partially en- 
tangled pure state into a maximally entangled state with the 
same dimension providing that they are close enough to each 
other. 

Theorem 6.1: Let \tp) be a state in S°(\ip)) and let |$+) = 
Si=i 7)t ^ e a ^ x & maximally entangled state. Then 
there exists a positive number 6 such that for any k x k state 
|x) satisfying |||x) — < 8, the transformation of |^>)<8>|x) 
to \p) <g> \& + ) can be realized with certainty by LOCC, i.e., 

\i>)®\x) -» 

Proof. This is a simple application of Lemma I2.3I Since 
S°(\ip)) is not empty and l u {\$ + )) = 1 > gu(\f)), it follows 
from Lemma |2"31 that 

ip (g) $ + < <p <g> <i> + . (41) 

An arbitrary but small enough perturbation on |$ + ) in the 
left hand side of Eq. (HTl i can still keep the relation '<'. 
Hence the existence of 6 is proven. □ 

The above theorem tells us that for any given \ip) <E S°(\ip)) 
and k > 1, we can find a partially entangled pure state 
|x) satisfying \ip) eg) |x) — > \<p) ® \® + )- It is obvious that 
|x) depends not only on \p) and k, but also \ip). At first 
glance, this seems to be contradicting our result about partial 
entanglement recovery, which states the auxiliary state |x) for 
partial entanglement recovery only depends on the target state 
\<p) and the presence of equalities in the majorization ip -< p. 
The key point is when we consider whether |x) can be used to 
do partial entanglement recovery for a transformation with the 
target state \<p), the resulting state \ui) is not specified; while 
the resulting state here is given and is maximally entangled. By 
Nielsen's theorem, |x) should be determined by the relation 
ip ® X V ® which obviously depends on the source 
state, the target state and k. 

Theorem 16. 1 1 confirms the existence of the partially entan- 
gled state |x). But it cannot yield a complete characterization 
of |x). To obtain such a characterization, we need to apply 
Nielsen's theorem and solve the corresponding majorization 
relation directly. To illustrate this procedure better, let us 
examine a simple case where \ip) and \<p) are both 2x2- 
dimensional. In particular, the following example deals with 
the case of k — 2. 

Example 6.1: Let \tp), \<p), and |<£> + ) be three 2x2 states 
with ip = (a, 1 — a), tp = (6, 1 — b), and <1> + = (|, ~), where 
| < a < b < 1. We are going to find a 2 x 2 partially 
entangled state |x) such that the transformation of \ip) to \<p) 
can concentrate \\) into the maximally entangled state |$ + ). 

Suppose that |x) is of the form \ — (j>A ~ p)< where 
h < p < 1. By Nielsen's theorem, we only need |x) to satisfy 

i>® x ■< * + - 

Notice that p® $ + has only two distinct components \b and 
1(1 — b). By the definition of majorization, the above equation 
holds if and only if 

ap < —b 
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and 



Hence 



(l-o)(l-p) >-(!-&). 



1 . b l-2a + 6, 

_ <p<min{ _ ___}. 

Note that | < a < 6 < 1, so the above equation can be 
simplified into 

\<P<^~, (42) 
I la 

which is exactly the result obtained in [10]. □ 

More generally, suppose that the k x k auxiliary state \x) is 
of the form \ = (71 , • • • , Ik) ■ Then to obtain a fc x k maximally 
entangled state |<1> + ) from the above transformation of \ip) to 
I ip), it suffices to have ip eg) x ~< f ® , which is equivalent 
to 



and 



Thus 



71 a < 



7 fe (l - a) > 



1-6 



1-6 6 

Tjz r < Ik < 71 < tt- 

fc(l — a) ka 



(43) 



If k = 2, we can show that Eq. d43l can be reduced to Eq. 
dHJ. 

One can similarly consider the general case where both k 
and n are arbitrary positive integers. 

B. Partial entanglement recovery and quantum catalysis 

In the above discussions, we always assume that the source 
state \ip) is comparable to the target state \<p), i.e., the 
transformation of \tp) to \<p) can be realized with certainty 
under LOCC. How about the case where \ip) and \tp) are 
not comparable? The general answer to this question remains 
unknown. 

In [14] a special case where the transformation of \ip) to 
\<p) has a catalyst state \c) such that \ip) ® \c) — > \<p) (g> \c) 
[18], i.e., the transformation of \ip) to \<p) can be realized 
under ELOCC, was examined carefully. It was shown that 
the problem of doing partial entanglement recovery for the 
transformation of \tp) to \<p) with ip -fi ip may be reduced to 
the problem of finding a catalyst state \c) and then seeking a 
suitable auxiliary state \x) to do partial entanglement recovery 
for the new transformation of \ip) ® \c) to \<p) ® |c) such that 
ip®c< (p(g)c. For this purpose, in [14] an algorithm of time 
complexity 0((nk)l) was proposed to find a k x k catalyst |c) 
for a transformation of \ip) to \<p) in which the source state 
and the target state are both n x n-dimensional. 

However, the above algorithm is intractable since it is 
of exponential time complexity. In [22] a polynomial time 
algorithm of n for fixed k was given. With the aid of this 
efficient algorithm, one can quickly determine whether an 
nxn incomparable pair has akxk catalyst. Then by the results 
obtained in the present paper, such as Theorems 13.21 14.11 or 
14.21 a state \x) that can do partial entanglement recovery for 
the transformation of \ip) to \<p) can be explicitly constructed. 



Therefore, if the transformation of \ip) to \ip) can be 
realized with certainty under ELOCC, then we can find an 
auxiliary state \\) to do partial entanglement recovery for this 
transformation. 



C. Partial entanglement recovery and mutual catalysis 

In [23], an interesting phenomenon named mutual catalysis 
was demonstrated. If \ip) \p) and |a) \J3) but \ip) ® 
\ a ) ~ ¥ \v) ® we say that \tp) and \a) can be mutually 
catalyzed by each other. The trivial case such that \ip) — * \/3) 
and \a) — > \ip) is not necessary to consider. With the help 
of the results obtained in previous sections, one can easily 
construct many non-trivial instances with the mutual catalysis 
effect. First, let us reexamine an example from [23]. 

Example 6.2: Let \tp), \ip), \a), and \(3) be four states with 
ip = (0.33,0.32,0.3,0.05), ip = (0.6,0.2,0.14,0.06), a = 
(0.6,0.3,0.1,0), and (3 = (0.46,0.46,0.08,0). It is easy to 
see that both the transformations of \tp) to \tp) and of \a) to 
\(3) cannot happen with certainty even under ELOCC. But we 
do have ® |a) — > \<p) ® \j3) in a non-trivial way. This is 
just the effect of mutual catalysis. 

From another point of view, this example can be 
treated as a special instance of partial entanglement 
recovery. To see this, let us relabel the above four states 
as follows: X = (0.6,0.3,0.1,0), u> = (0.6,0.2,0.14,0.06), 
ip = (0.33, 0.32, 0.3, 0.05), and ip = (0.46, 0.46, 0.08, 0). It is 
obvious that ip <\ ip. Noticing that l u (\ui)) > g u (\f)) = 0, we 
have that ip ® us <\ p ® oj by Lemma l2~3l A small perturbation 
on u) will generate \ = w + (0, 0.1, -0.04, -0.06). Note that 
E(\x)) = 1.2955 < E(\u)) = 1.5472. So the entropy of 
entanglement of \x) is enhanced. □ 

The above example suggests a connection between partial 
entanglement recovery and mutual catalysis. More generally, 
any pairs {\ip),\x}} and {|</?),|w)} such that \ip) <g> \x) -> 
\<p) ® |x) -» |x) \<p), and \ip) ^ \p) (g \uj) can be 
treated as nontrivial instances of mutual catalysis. These pairs 
can be easily obtained with the aid of lemma 12. 3 1 Furthermore, 
one can choose the state \uj) satisfying u) <x but X -fi u - We 
omit the construction details. 



D. Multiple-copy is essential for partial entanglement recov- 
ery 

Multiple-copy entanglement transformation is another in- 
teresting topic in quantum entanglement theory. Let us review 
this concept briefly. In [15], it was demonstrated that some- 
times multiple copies of a source state may be transformed 
into the same number of copies of a target state although the 
transformation cannot happen for a single copy. That is, for 
some states \ip) and \tp), although the transformation of \ip) 
to \<p) cannot be realized with certainty by LOCC, there may 
exist m > 1 such that the transformation of \ip)® m to \ip)® m 
can be achieved deterministically. This kind of transformation 
that uses multiple copies of a source state and then transforms 
all of them together into the same number of copies of a 
target state is intuitively called 'multiple-copy entanglement 
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transformation', or MLOCC for short. See [15], [19], [20], 
and [21] for more about MLOCC. 

It may be of interest to study the relations between partial 
entanglement recovery and multiple-copy entanglement trans- 
formation. To our surprise, entanglement lost in a multiple- 
copy entanglement transformation can be recovered more 
easily than that in a single-copy transformation when the 
auxiliary state is specified. To demonstrate this point, we need 
the following theorem as a useful tool. 

Theorem 6.2: Let \<p) and |x) be two partially entangled 
states. If \x) has at least two distinct nonzero Schmidt coef- 
ficients, then there exists a positive integer ko such that for 
any k > ko and \ijj) £ S°{\ip)), entanglement lost in the 
transformation of \ip)® k to \(p}® k can be partially recovered 

by|x>- 

The most interesting part of the above theorem is that the 
choice of k only depends on |%) and \ip). 

Proof. First, applying Lemma 1 in [20] yields that rp < <p 
implies ^® fe < ip® k for any k > 1. Second, notice that 
9u{\<p)® k ) = fl£(|<P» and g u (\<p)) < 1. By the assumption 
on \x), we have < L u (\x)) < 1. Thus there exists ko > 1 
such that L u (\x)) > 9u(\f)) f° r anv k > ko- Therefore, by 
Theorem 13.21 we deduce that |x) can be used to do partial 
entanglement recovery for the transformation of \ip)® k to 
\ip)® k such that k > k and \ij>) £ S°{\tp}). □ 

Let us take now two states \<p) and \x) such that < 
-Mix)) < 9u(\v)) < 1- By Theorem |3T2l |x) cannot do 
partial entanglement recovery for any transformation with the 
target \<p) since L u (\x)) < g u {\(p)). On the other hand, 
it is easy to see that \<p) and \x) satisfy the assumptions 
of Theorem 16.21 Hence, there exists ko such that for any 
k > ko and \ip) £ S°(\(p)), entanglement lost in the fc-copy 
transformation, i.e., the transformation of \ip)® k to \<p)® k , can 
be partially recovered by |%). 

We give an intuitive explanation for the above theorem. If 
the auxiliary state \x) cannot do partial entanglement recovery 
for the transformations with the target state | ip) , then the target 
state is too uniform, and it is too entangled relative to |x)- 
So for any state that can be transformed into \cp) under 
LOCC, the extra entanglement left (except the necessary part 
to finish the transformation of \tp) to \(p)) is not enough to be 
transferred into |x). But if multiple copies of the source state 
are provided, the extra entanglement will accumulate. Such 
extra entanglement can be transferred into the state \x) when 
it exceeds a threshold. 

It is also interesting to investigate the partial entanglement 
recovery power when multiple copies of the auxiliary state 
\x) are available. We restrict ourselves to the special case that 
|x) is 2 x 2-dimensional. A surprising result appears as the 
following: 

Theorem 6.3: Let |x) be a 2 x 2-dimensional partially 
entangled state and \ip) be any partially entangled state. If 
L u (\x)) < 9v(\<p))> m en for any k > 1, \x)' g " k cannot 
do partial entanglement recovery for any transformation with 
target state \ip). 

Intuitively, if the auxiliary state |x) is 2 x 2-dimensional, 
then more copies of |x) do not provide any extra power of 



partial entanglement recovery if L u (\x}) < gu(\<p))- This is 
very reasonable since the key point of such a recovery is the 
difference between the entanglement resource of the source 
state and that of the target state, which keeps invariant during 
the process of increasing the number of copies of the auxiliary 
state. 

Proof. By (3) of Lemma 12.21 we have that 
Lu(\x)® k ) = -Mix)) since \ X ) is a 2 x 2 state. So by 
the assumption L u (\ip)) < g u (\ip)), it follows from Theorem 
13.21 that |x)® fe cannot be used to do partial entanglement 
recovery for any transformation with the target state \<p). □ 

In the case when L u (\x)) = 9u{\f)), however, the partial 
entanglement recovery capability of \x)® k may be strictly 
more powerful than that of |x) for suitably large k. That 
is, |x)® fe can do partial entanglement recovery for some 
transformation with the target state \<p) while |x) cannot. See 
the following example. 

Example 6.3: Let \x) be a state with x = (p, 1 — p), where 
2 < p < 1. Then by Theorem 13.21 we know that |x) cannot 
do partial entanglement recovery for any transformation with 
target state \<p) such that ip = (p,p,p,p, 1— p, 1 — p)/(2 + 2p). 

However, by Theorem 13.21 again, |x)® 2 with 
X® 2 = (p 2 ,p(l — p),p(l — p), (1 — p) 2 ) can do partial 
entanglement recovery for any transformation of to \<p) 
such that \ip) is in S°(\(p)). □ 

A more general result in this special case is: if \ip) has 
only two distinct nonzero Schmidt coefficients, then for a 
sufficiently large k, \x)® k can always do partial entanglement 
recovery for any transformation of \ip) to \ip) with \ip) £ 
S°( | </?)); otherwise such a recovery is impossible for arbitrarily 
large k. 

VII. Conclusion 
To summarize, we obtain a complete characterization of 
an auxiliary bipartite entangled state |x) that can do partial 
entanglement recovery for the transformation of to \<p) 
where if) is strictly majorized by ip. It is interesting that the 
choice of the auxiliary state can only depend on the target 
state \ip) and the presence of the equalities in the majorization 
relation i/j -< (p. We further propose two sufficient conditions 
for |x) that can be used to do partial entanglement recovery 
for a class of transformations of \ip) to \<p) with t/j -< (p. We 
also study the feasibility of partial entanglement recovery from 
the algorithmic viewpoint. A polynomial algorithm of time 
complexity 0(n 2 k 4 ) is presented for deciding the possibility 
of partial entanglement recovery. As applications, we establish 
some interesting connections of partial entanglement recovery 
to the generation of maximally entangled states, quantum 
catalysis, mutual catalysis, and multiple-copy entanglement 
transformation. We hope the results presented here may help 
us to manipulate quantum entanglement more economically. 

Appendix A: Proof of Lemma 12.31 
Take \ip) £ S°(\<p)). Assume ijA = (ori, a2, ■ ■ ■ , a n ), = 
{01,02, ■ ■ -,0n)< and x l = ill, 72, ■ ■ ■ ,1k)- If k = 1 then the 
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result follows trivially. In what follows we assume k > 1. 

First we prove that if l u (\x)) > 9u(\<f)) then \if>) ® |x) is 
in S°(\(p) <g> |x))- In other words, we shall prove 



Appendix B: Proof of Theorem|3.2| 



for any 1 < I < nk. 
We rewrite 



(44) 



(45) 



where < U < n and Ei=i U = Easily see that 7^ < 7^93 
for all 1 < i < k. So we merely need to consider the following 
two cases: 

Case 1: There exists 1 < s < k such that < l s < n. In 
this case, e,i B (lsi>) < ei s ilsf) holds. Then Eq. (l44l follows 
from 

ei(ip®x) = E;=i e h (TiVO 

< Ei=l e h ilif) 

< e t {(p®x), 

where the equality is by Eq. d45l l. and the second inequality 
is by the definition of ei(tp £g> x). 

Case 2: For any 1 < i < k, U 6 {0, n.}. Let /i be the 
maximal index satisfying i/, = n. Then 1 < h < k; otherwise 
h = k implies I — nk, which contradicts the assumption I < 
nk. Noticing l u (\x)) > 9u(\<P)), we have •Jh+i/lh > Pn/Pi, 
or 

IhPn < 7h+i0v (46) 

By the definition of e n h (tp <8> x) an d me assumption on h, we 
further have 

e n h(i> ® x) = Ei=i e n{li^) 

= EiJl e„(7iy) + E?=l IhPi + IhPn, 

where the second equality is due to e n (ip) = e n (ip). 
Substituting Eq. (l46l l into the above equation, we have 

e n h(ip ® X) < EiJi 1 e «(7^) + Ei^i 1 IhPi + Ih+iPi 

where the second inequality is by the definition of e n h{<P®x)- 
Therefore Eq. (144") holds for any 1 < I < nk. By the 

definition, we have ® |x) is in S°(\ip) ® |x))- (Note that 

e n k{ip ® x) = e nk(f <S> x) holds for any x whenever ^ < tp). 
Conversely, suppose ® |x) is in 5 , °(|(/?) ® |x)), while 

there exists some 1 < h < k, such that 

Ih ~ Pi ' 
IhPn > Jh+lPl- 

e n h(<P ® X) = Ei=i e n (7^) 
= EiU ^(^VO 

which contradicts the assumption that ei(ip eg) x) < e i{f ® X) 
for any 1 < I < nk. 

With that we complete the proof of Lemma 12.31 



or equivalently, 



Then we have 



To prove Theorem l3.2l we need the following three auxiliary 
facts about majorization. 

Lemma 7.1: Let y E TZ n with compact form y+ — 
(y'® kl , ■ ■ ■ , y' s ® ks ), and let x e TL n satisfying x -< y but 
y -/< x. Then there exists z <G lZ n such that x -< z -< y 
and 



r 'fflfci /ffifci-l / 

^ = [j/i \ • ■ • , Vi , V\ ~ e, 



1 . i®kj-l 



,y 



for some 1 < i < j < s and e > 0. 

Proof. This is a direct consequence of B.l. Lemma in [7] (page 
21). □ 

Lemma 7.2: If x' -< y' and x" -< y" , then (x',x") -< 

(y',y")- 

Proof. See part (i) of A.7. Lemma in [7] (page 121). □ 

Lemma 7.3: Let — (x'^,x"^) and y^ — {y 1 ^, y"^). If 
x -< y and x' -< y' (or x" -< y"), then x" -< y" (resp. 
x' -< y'). 

Proof. Suppose x',y' € K m and x",y" S 1Z n . By the 
assumption, we have 



and 



ei(x") = e m+ i(x) - e m (x') 



ei(y") = e m +i(y) - e m (y') 



for any 1 < I < n. 

Noticing x' -< y' and x < y, we also have 



and 



e m (x') = e m (y') 



ei(x) < ei(y) 



(47) 
(48) 

(49) 
(50) 



for any 1 < I < m + n and the inequality is an equality if 
l = m + n. Thus Eqs. g7)- (|50) give 

e l {x")<e l {y") 

for any 1 < I < n, with the equality holding when I = n. 
That means x" -< y". If x -< y and x" -< y", we can prove 
x' -< y' similarly. □ 



Now we proceed to the proof of Theorem 13. 21 We first deal 
with the case where all the nonzero Schmidt coefficients of 
|x) are identical, i.e., L u (\x)) = or L u {\\)) = 1. 

Assume x 1 = ((^)® a , e/c_a ). If a = k, then |x) is a k x fc 
maximally entangled state. There cannot be another k x k state 
|<a>) that is more entangled than |x). So partial entanglement 
recovery is not possible in this case. 

Now suppose 1 < a < k. If |x) can do partial entanglement 
recovery for some transformation of \ip) to \tp) such that \ip) 
is in S(\(p)), then the state \ui) such that rp ® x ^ <^ (8> ^ and 
w -< x> but x^ 7^ should have at least a+1 nonzero Schmidt 
coefficients. By the property of majorization, the number of 
nonzero Schmidt coefficients of ip ® x i s not l ess man mat °f 
ip ® to, i.e., 

n"a > n'(a+ 1), 
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where n" denotes the number of nonzero Schmidt coefficients 
of \ip). Obviously, n" < n, thus we have 

na>n'(a + l). (51) 

Conversely, if Eq. (IBTl i holds, we will show that |x) can be 
used to do partial entanglement recovery for any transforma- 
tion of \ip) to \<p) with \if>) £ S°(\ip)). Let us take 

X ( c ) = ((l)e-i,i_ e>e , ffi*— i), 
a a 

where e is a suitably small positive number. We have that 

X(e) -< x and x l + X X ( e ) 

for any < e < ^. Noticing \ip) € S°(\tp)),by Lemma |231 
we have 

(^®(-)®°)<(y®(-) ffl °)- 
a a 

Combining the above equation with Eq. (IBTl i gives 

® V°>0))i o , 

a a 

where the notation (x)* denotes the segment (x\, . . . ,x\) of 
ar. Since a sufficiently small perturbation on the right-hand 
side of the above equation cannot change the relation <, we 
have that 

(iP ® (-) 0a ) < (<P ® ((-) 0Q - 1 , - - e, e))^ (52) 
a a a 

for small enough positive number e. 

Appending suitable number of zeros on the both sides of 
Eq. ( |52l gives 

V> ® x ^ v ® x( £ )i 

which proves that |x) can do partial entanglement recovery 
for the transformation of \ip) to \<p). 

Now we turn to the general case that |x) has at least two 
nonzero distinct Schmidt coefficients, i.e., < L u (\x)) < 1. 
We shall consider the following three cases: (i) L u (\x)) > 
g u {W))\ (ii) -Ml*}) < 9u(\f)); and (iii) L„(|x» = 9u(\<f))- 

First, we deal with case (i). We shall prove that if L u (\\)) > 
g u ( | ip)), then |x) can do partial entanglement recovery for the 
transformation of \ip) to \ip) such that \ip) is in S°(\ip)). 

Suppose x l = (7?*V-->7m* m )- Then there exists 1 < 
i < m such that 

gM) < L u (\x)) = 2*±1 <i. (53) 

7» 

Let us denote (7?* 1 , • ,7®** -1 ). (7i,7i+i). and 
(7® fe r l_1 ,...,7® fcm ) by 7', 7", and 7"', respectively. If 
ki — 1 or fcj+i = 1, we simply omit the meaningless terms 
7' or 7"', respectively. For any \ip) G S l °(|(^)), we have 

V>® 7' -<! <p ® 7', (54) 

V> ® 7" < <p ® 7", (55) 

and 

^ ® 7"' -< V ® 7"', (56) 



where Eq. (|55] l comes from Eq. ( 1531 and Lemma [23] So for 
a sufficiently small positive number e, we have 

ip ® 7" < ip ® 7"(e), (57) 

where 7 "(e) = (7* - e, 7^+1 + e). By Eqs. dSU, (|56]) and 071 
applying Lemma 17.21 gives 

"0 ® X ~< ® W) 

where x 1 = (7', 7", 7") and u l = (7', 7"(e), 7"'). It is 
obvious that u -< x but x^ 7^ w^. 

Second we deal with case (ii). Suppose that L„(|x)) < 
<7„( I ip)). We shall prove that |x) cannot do partial entangle- 
ment recovery for any transformation — > \ip) such that \ip) 
is in S(\(f)). 

By contradiction, suppose that there exists a state \u>) such 

that ip®x~4'P® u > w ~4X an d X^ 7^ ^ ■ F° r an Y state lx') 
such that u> -< x' ~< X we have 

(58) 

where we have used the assumption ?/> ® x ^ V ® w - 
By Lemma [7TT1 x' can be chosen as 

x' i = (7',7"W,7'"), (59) 

where 7' = (7® fcl , . . . , 7 ® fc *- 1 ), 7 "(e) = ( 7i - e, . . . , 7i + e), 
rf" = ( 7 ® fc ^ 1 ,..., 7 © fe ™), 1 < i < j < m, and e is an 
arbitrarily positive but small enough real number. In particular, 

X i = (7',7",7"'), (60) 

where 7" = 7 "(0) = (7,, . . . , 7,). However, we shall prove 
that such two indices i and j cannot exist, and thus complete 
the proof of this case. 

For simplicity, let n = dim(<^). By the assumption 
Mix)) < 9u{\v)), it follows that 

j s f3 n > 7 s +i/3i, for any 1 < s < m - 1, (61) 

where j3\ and /3„ are the greatest and the least components 
of ip, respectively. Notice that ip -< (p. Applying part (6) of 
Lemma Ogives g u (\ip)) > g u {\<f)), thus L u {\x)) < 3«(|V>))> 
or more explicitly, 

7 s a„ > 7s+ iai, for any 1 < s < m — 1, (62) 

where a% and a n are the greatest and the least components of 
ip, respectively. 

Eqs. d60b and d62b imply 

(^(gjx) 4 = ((0«)7') i ,(V'«)7") i ,(V'«)7'") i )- (63) 
Eqs. j59l ), d6TT l. and e > imply 

(ip^x') 1 = ((<p®y) l ,(<P®l"(t)) l ,(<P®l / ") 1 )- (64) 
Applying Lemma |7~3l to Eqs. d63l and ( |64] i yields 

0®7" -< V® 7 "(e), (65) 

where we have used the assumption ip -< y> and Eq. (|58l . 

According to Eq. doTT l. we can take a sufficiently small 
positive number e such that 

(7i - e)/3„ > (7^+1 + e)/3i. 
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Then 

e n (f ® 7"( e )) = {li ~ e ) e n(<p) 
< e n (ip®y), 

which contradicts Eq. d65b . 

Finally, we deal with case (iii), i.e., L u (\x)) — 9u(\<p))- 
This case is much more complicated than the previous two 
cases. It is in fact the most non-trivial part of Theorem 13.21 
Since this case is of considerable interest, we will present a 
detailed proof for it. To keep the proof as readable as possible, 
the lengthy proof is divided into two easier lemmas. It is worth 
noting that both lemmas are interesting in their own right. 

The first lemma shows that an auxiliary state can do partial 
entanglement recovery for a specific transformation if and only 
if some of its segments can do partial entanglement recovery 
for the same transformation. 

Lemma 7.4: Let \ifj) and \<p) be two states such that if) -< 
if, and let |x) be a partially entangled state with compact 

form x l = (7? fc V ■ • ,1m ) for some m > 1. If L u (\x)) = 
g u (\ip)), then the following two statements are equivalent: 

(i) |x) can do partial entanglement recovery for the trans- 
formation of to | ip); 

(ii) There exists an index i such that the unnormalized state 
\x') with x'^ — {lf ki i 7i+i +1 ) can do partial entanglement 
recovery for the transformation of to \<p), where 1 < i < m 
and 2^1 = g u (\tp)). 

Proof. The essential part of the lemma is (i) =>■ (ii). 
Suppose that \x) can do partial entanglement recovery for the 
transformation of \ip) to \ip). That is, there exists a state \u) 
satisfying ?/><g)x^</? ( 8)w, w^x> and X 1 7^ w ■ Moreover, by 
Lemma 17721 we can assume that \u) is of the following form: 

^ = (y,y'( e ),y"), 

where y = ( 7 ? fcl , . . . , 7^ ). 7"(e) = (it^^li - 

e,- ...7i + e,7f Y" = • ..tS*"). 1 < * < 

j < m, and e > 0. To make (l"(e)Y = '/'(e) hold, we 
have assumed that e satisfies 7 P — e > 7 P +i + e for any 
1 < p < m — 1. We also have 

X i = ( 7 ',7"(0),7'"), 

where 7 "(0) = ( 7 r%..., 7 f0- 

By the assumptions L u (\x)) = 9u(\*p)) and ijj -< ip, we 
have 

(V-^x) 4 = ((V'®7') i ,(V'«)7"(0)) i ,(V'®7"') i )- (66) 
and 

[if <g> uj) 1 = ((ip <g) 7 ') X , (y> ® 7 "(e)) X , (y> <8> 7 "') X )- (67) 
Applying Lemma 17731 to Eqs. d66i l and (|67T > yields 

^® 7 "(0) -< Lp®i'(e), (68) 
where we have used the assumptions that tp ® x -< ¥> <£> u> and 

Therefore, for the simplicity of notations and without any 
loss of generality, we can assume that i = 1 and j = m. 
More directly, we can write 7 "(0) and 7 "(e) as x and u, 
respectively. This, of course, will not cause any confusion. 



We shall prove m = 2 and ^ — ffud'f 5 )) to complete the 
proof of the lemma. 

For simplicity, we assume that n = dim(tp) in the rest of 
proof. 

First, we prove 1 < m < 3. By contradiction, suppose that 
m > 3. Let us decompose 

x = (x'(Q),x"(o)) 

and 

u ) = (x'(e), X "(e)), 

where 

X (e) = (7i »7i - e,72 ) 

and 

X"( e ) = (73 efc3 ,---,7 m + e,7^ 1 ). 
Again, L u (\x)) = and V ^ <^ give 

73/72 < Su(lv)) and 73/72 < 5u(l^))- 
That immediately yields 

(-0 (gi x) X = O x'(o)) x , (V> ® x"(o)) x ) 

and 

(if. <g> w)l = (foj <g> x'(e)) X , fo> ® x"(e)) X )- 

So, 

e (fel+fe2 ) n (V> ® x) = e(fe 1 +fe 2 )«(V' ® X'(0)) 

= fcl7l + fc 27 2 

and 

e(fc 1 +fe 2 )«(^ ® w) = e (fcl+fe2) „((p <g> x'(e)) 
= fci7i + /c 27 2 ~ e, 

thus 

e(/c 1 +fe 2 )n(V' ® x) > e(fc 1+ A a )nOP ® w ) 

for any small e > 0. This contradicts the assumption 7/1 eg) x ~< 
tp®u). Hence 1 < m < 3 

Second, we prove that for any 1 < i < m — 1, it holds 
that 7 ;+i/ 7 i = P n /I3i, where /?i and /3„ are the greatest and 
the least components of ip, respectively. By contradiction, we 
need to consider two cases: (1) m = 2 and (2) m = 3. 

(1) m = 2. Suppose that 72/71 < (3 n / Pi- Let us choose a 
suitably small positive number e such that 

(72 + e)/?i < (7i -e)P n - 
A routine calculation shows that 

efcmOOx) = e / c 1 „0®7i e ' £l ) 
= fci7i 

and 

e kin (<p ®uj) = e (fcl _ 1)n ((/3 ® 7i efcl_1 ) + e„((7i - e)v) 
= fei 7 i - e, 

which yields 

e/fcmC^ ® X) > efc in (^ ® cj) 
for any small e > 0. That again contradicts t/j ® x ^ ® w - 



21 



(2) TTi — 3. Suppose that 72/71 < 8 n / j3 x or 73/72 < 
Pn/Pi- We only consider the case where 73/72 < fln/fli, 
and the left case is similar to case (i). Choose a suitably small 
positive number e such that 

(73 + e)A. < 72/?™ • 
Then a simple analysis shows 

e ( fe 1 +fc 2 )n(V' ®X) = efcmO <8 7? fcl ) + e fc2 n(V> ® 7® fc2 ) 
= fci7i + fc 2 72 

and 

e(fc 1 +fc 2 )n(<^ ® w) = e( kl _i)„(¥) ® 7? fel_1 ) 

+e„((7i -e)^)+e fe2 „(<^®7® fc2 ) 
= fei7i + A: 2 72 - e, 

which yields 

e(fc 1+ fc 2 )™W> ® X) > e( fel+fc2 )„0 (g) w). 

That is a contradiction with ?/> ® x -< f S3 u;. 

Third, we prove that m = 2. By contradiction, we shall 
show that if m = 3 then 

V 1 ® x ~< v ® w 

cannot hold for any small enough positive number e, where 

X 4 = (7i ^ 1)7 e fc2)7 e fe3)) wi = (7 e^-i )7l _ ^efc^ + 

e,lf k3 1 ). 72/71 = 73/72 = 

To be specific, let ^ = {df mi , . . . , (3® mh ) for some ft > 2. 
Obviously, n = dim(^) = J^ t=1 Tit- Choose e such that 

(7i - e)A > 7iA+i 

and 

(73 + < 73A 

for any 1 < i < ft — 1. 

In addition to the above constraints, we also need e satisfy- 
ing 

(71 - e-)lh > (73 +e)A- 
A direct calculation gives 

(V'®x) i = (^(o),^'(o)) 

and 

1 - r,„'iv> 



te® W )* = fe/(e),p"(e)) 



where 



and 



= [ 7l /3® (fcl - 1)m ' l+ ' £2mi , (71 - e)/?rN 

...,(73+ e )/3r ft ,73/3r 3 " l)mh ] ; 

where we have used 71/^ = 72/81. 
It can be readily verified that 

efc ltt O®x) = efc 1 (n-m h )(¥' / (0)) + e Aimh (vj"(0)) 
= fei7i 



and 

e kin (tp® u) = e kl{n _ mh) (<p' (e)) + e kimh (<p"(e)) 

= (fci7i - - rnhPh) + e k imh {tp"{e)). 

To calculate e^ mh ((f"(e)), we need to consider the follow- 
ing two cases: 

(a) k2Tn± > rrih- Then 

= kirrihjiPh, 

thus 

e kin {ip®uj) = fei7i - e(l - m h Ph) 
< e kin (%l>®x) 

providing e > 0. 

(b) fc 2 m 1 < m^. Then 

e fcim >"(e)) = e h ( 7 i/3f + e Ia (( 7 i - e)/?f 3 ) 
= hmhJiPh - e(m h - k 2 m 1 )(3 hl 

where l\ = (fci — l)irih + k^mi and Z2 = mh ~ fcjjni. Thus 

e fel „(^(g)w) = fci7i - e(l - femtft) 

< ^171 - e(l - mji/3/,) 

< e fcl „(-0 ® X) 

providing e > 0. In the above two cases we have used 7i/3/i = 
72/81 to simplify the calculations. 

Both the above two cases contradict tp <g> x ~< ¥> ® w. Thus 
m = 3 is impossible. 

With that we complete the proof of Lemma 17.41 □ 

By Lemma l7~4l under the condition L u (\x)) = 9u{\<f)), we 
only need to consider a special form of |%). More precisely, x 
has only two distinct components. The following lemma will 
just handle such a special form of \x). 

Lemma 7.5: Let \x) be a partially entangled state with 
compact form x l = {Pf kl , Pf 2 ) for some Pi > 02 > 0, and 
let \ip) be another state satisfying g u (\ip)) = L u (\x)). Then \x) 
can do partial entanglement recovery for the transformation of 
\ip) to \ f) such that \ip) is in S°(\ip)) if and only if 



(69) 



where x' is a segment of X' 1 with two distinct components, C 
is a normalization factor, and m > 1. 

Moreover, if |x) and |<p) don't satisfy Eq. j69l , then \x) 
cannot do partial entanglement recovery for any transformation 
of to \<p) such that \ip) is in S(\ip)). 

Proof. We first prove that if |%) can do partial entanglement 
recovery for some transformation of to \ip) with ip -< 99, 
i.e., there exists a state |u/) satisfying V , ®X~^¥'®' A '« W ~^X> 
and 7^ w^, then |x) and should satisfy Eq. 

Suppose that \tp) and \u>) have compact forms 

<P l ={Pf m \ 



and 



,1 - 



(7i ,7i -e,72 + e,7? )■ 



We shall prove that if ip x ^ 'Z 3 ® w anc l w ^ X f° r an Y 
sufficiently small positive number e, then 



1 TTii 

ft = 2 and — < < fci. 

k 2 m 2 



(70) 
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Or more compactly, <p has the form as in Eq. d69| i. 
The condition L u (\x)) = g u {\<p)) is equivalent to 

72/?i=7i^- (71) 
For any 1 < i < h — 1, choose e such that 

(71 - e)Pi > 7 i/3i+i, (72) 

and 

(72 + e)/3 i+1 < 7 2 ft. (73) 

In addition to the above conditions, we also choose e 
satisfying 

(7i-e)j8h-i >(72+e)A, (74) 



and 



(71 - > (72 + e)p2 



(75) 



By the condition L„(|x)) = gu{\f)) and V> -< <A it is easy 
to verify that 

(V, ® X ) 4 = ((V ® 7 f fcl ) i , W ® 7 2 ffifc2 ) i )- (76) 



Take 



dim^) = '< 



i=l 



Then by Eq. 

ek in (ip <S> x) = fci7i- 
By Eqs. (TTTTt — (T73T>, a careful analysis gives 

where 

*/ = [ 7l ^( fe - 1 ),( 7l -e)/3r i , 



(77) 



(78) 



and 



So 



•>7lPfc-i »(7i-e)P h _i ] 



( 71 -e)/3rN---,7 2 /3r h(fc2 ^ 1) ]- 



(80) 



e kin (<f <E) uj) = e kl{n _„ lh) (ip') +e kimh ((f") 

= (fci7i - e)(l - rrihPh) + ei 1% (/). 

(79) 

We need to consider the following four cases according to 
the values of e kinih (ip"): 

Case (a): mi > k\mh- Then it is obvious that 

e feim >") ^e kimh (( l2 +e)(3f k ^) 
= kimh{j2 + 
Combining Eqs. d79l with (1801 . we have 

efc irl (<p (g) w) = fci7i + e{k 1 m h l3i + m h f3 h - 1), 

where we have used the relation 71 p\ — I2P1 to simplify the 
calculations. 

Since V>®X ~< f^ui, it follows that e kin ((p®w) > e kin (ip® 
X), i-e., 

fci7x + (.{k-imufi-y + mftpVj - 1) > fci7i- 



Or equivalently, 

fcim h( Si + m/A > 1. (81) 
However, by mi > fcim/i and 53i=i = 1, it follows that 

hmnPi + m h (3 h < mi/3i + m h (3 h < 1, 

which contradicts Eq. d8 1 1 > . 

Case (b): < mi < k\mh. It is easy to calculate that 

e kimh (v") = emi ((72+e)/?r i ) 

+ e kimh - mi {liPT imh ~ mi ) ( 82 ) 
= emi/3i + kimhjxPh. 

By Eqs. d79j and ([82j, it follows that 

e kin {v> ®u) = fci7i + e(mi/?i + rn/j/3;, - 1). 

Since ?/><8>x -< y(8)u;, it follows that e kin (<p®uo) > e kl n(ip® 
x)> i-e-, 

fci 7 i + e(mi/?i + mfoflii - 1) > fci7i. 
Or equivalently, 

roi/?i + mjA > 1. (83) 

It is easy to verify that Eq. d83l holds if and only if h = 2. 

Case (c): mi < m/j < k 2 m\. Similar to Case (b), ip ® x ^ 
(p ® ui holds for any small enough positive e if and only if 
/i = 2. 

Case (d): m^ > k 2 m\. Similar to Case (a), this also causes 
a contradiction. 

Summarizing the above four cases, we obtain that \(p) 
should satisfy Eq. (l70l . which is equivalent to Eq. d69l . 

Now we turn to prove that the condition in Eq. d69l is also 
sufficient for partial entanglement recovery. Suppose that \<p) 
and \ui) are with compact forms 

- i = (7r i - 1 ,7i- e ,72 + e ,72 ffi ^ 1 ), 



and 



where 



P2 72 , 1 , mi , , 
— = — and — < < Ki. 

pi 71 k 2 m 2 



(84) 



Take \ip) £ S°(\ip)). We shall prove that for a sufficiently 
small positive number e, the transformation of |x) to 
\ip) (gi \u>) can be realized with certainty under LOCC. 

By the assumptions L u (\x)) = 3«(|v)) and tp -< <p, it is 
easy to verify that 

1^^ = ^/), (85) 
where ip' = (ip ® 7f fcl ) x and ?/>" = (V> ® 7® fc2 ) i - Similarly, 

(^(ax)^^,^'). (86) 

where = ((/?® 7 f fel ) J - and ip" = (tp®jf k2 ) 1 . 

By Eq. ( f84b . it holds that 71/^2 = 72/^1 • Hence we also have 



and 



^ = (7i/?i ffifcimi ,72pf fcim2 ) 



^ = (7i ^efe 2 m 1)72/3 ew 2) . 



(87) 
(88) 
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Similarly, 



where 



(89) 



and 



II \ r a®(kl — l)mi / \/3©mi 

V(e)= [7i/3i >(7i-e)Pi , f9(r) 

(72+e)/3f mi ; 72/3f fcim2 ^ mi ] 

= [ 7 i/3 2 ffiroifc2 - m2 ,(7i-e)^ ero2 , 

( 72+e )/?r 2 , 72/3 2 e(fc2 ~ i)m2 ]. 

Note that Eqs. < f90b and (|9TT > are well-defined since we have 
Eq. (|84| |. We have also assumed that e in Eqs. < f90b and (|9"TT > 
satisfies the following constraints: 

(71 - e)ft > ( 72 + e)/3i and (72 + e)(3 2 < (71 - e)&. (92) 



Since ?/><](/?, by Eqs. ( |85l > and ( I86I 1, applying Lemma 12.31 
gives 

?// < 93' and V" < ¥>"• (93) 



A careful observation caries out that (f'(e) and <p"(e) are 
obtained by adding perturbations on if' and ip", respectively. 
So we have 

ip' < y/(e) and V" < p"(e) (94) 

for a sufficiently small positive number e. 

Thus by Eqs. d85l l. (|89l , and (|94l , applying Lemma l7/2l gives 



V"8> X f ® UJ - 



(95) 



It is easy to see that u> -< x ar, d X^ 7^ w providing the 
positive number e small enough. In other words, |x) can do 
partial entanglement recovery for any transformation of \ip) to 
\<p) such that \ip) is in S°(\<p)). 
With that we complete the proof of Lemma 17.51 □ 



Appendix C: Proof of Theorem I4.1I 

For simplicity, we denote and D$ u> by I and D, 

respectively. We only need to show that for any 1 < I < 
dim(-0)dim(x), one of the following two cases holds: 

Case 1: ei(ip <g> x) < ei(ip ® x); or 

Case 2: ei(ip <g> x) = e;(yj ® %), but both sides are not 
related to x* (i € -D) an d both of them remain unchanged 
by an arbitrary but small enough perturbations on x*(* £ -D)- 
Here we should point out that x % is supposed to be kept as a 
constant for each i E D during the perturbations to guarantee 
that |x) is a valid quantum state. 

For this purpose, we rewrite ei(ip <8> x) as follows: 



m m 



ei{ip <8> x) = 2^ 1^ e k,M l ® X J ), 
i=i j=i 

where J7 ■ l^j = I and < < dim^ ® x')- 
It is easy to see that 

< ei(ip(g>x), 
where the first inequality follows from ip l 



(96) 



If one of these inequalities is strict, then Case 1 holds, and 
the proof is completed; otherwise we only need to prove that 
Case 2 holds. 

More precisely, we only need to show that if ei (ip ® %) = 
ei((f®x) m en for any 1 < i < m and j 6 D, ijj can only take 
two values: or dim('i/; 4 <g>x : ') ■ Notice that IUD = {1, . . . , to} 
and 7 n D = 0. It suffices to prove E {0, dim(V> 4 ® X 5 )} 
for two cases: (1) i E I, j E D, and (2) i E D, j E D. 

Let us consider the case when i E I and j E D first. By 
Eq. (fTOt . we have 



</? 1 <8> x J C v 5,7 ® x\ for all i E I and j E D. 



(97) 



That is, the values of the extreme components of p % (g) x J are 
strictly bounded by those of pi ® x 4 - Thus, we have 

< (^®x 4 )i 

and 

(^0X 3 )I > (^®X*)l 
Hence by the assumption that all inequalities in Eq. ( |96*l l hold 
with equalities and the definition of ei((p ® x)> together with 
the above two equations, we have 







and 



= dim(V> J O X 1 ) h,j = dim(V>* O x J )- 

So, in order to prove Z.y S {0, dim( , I (g)x : ' )} m the case of i £ 
7 and j 6 £>, we only need to show that Ijj E {0, dim(^ (g> 
X 1 )} for i E I and j 6 Z?. Or equivalently, to show ijj 6 
{0, dim(-0 l <S> X 3 )} f° r i ^ D and j 6 7. (Here we interchange 
the indices i and j for convenience.) 

So combining this with the case of i E D and j E D, the 
only thing left to be proven is that 

kj E {0, dim^t&x )} for all i G D and 1 < j < m. (98) 
By Eq. (fTTT i and Lemma |231 we have 

ijf ® x 1 <\f' 1 ® X J for all i E D and 1 < j < m. (99) 

If there exist s E D and 1 < t < m such that < l S: t < 
dim(-0 s ® X*): then by Eq. (l99l we have 

e; 8>t (V' s ®X t )<e^(^®X t )- (100) 
It follows that the first inequality in Eq. ( |96l > strictly holds, 
which contradicts e;(^ <8 x) = e i(^ ® x)- So Eq. ( |98l holds. 
With that we complete the proof of Theorem 14.11 



and the second one follows from the definition of ei(<p (8) x)- 
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